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Abstract 



One of the most fundamental questions we can ask about a given gauge theory is its phase 
diagram. In the standard model, we observe three fundamentally different types of behavior: QCD 
is in a confined phase at zero temperature, while the electroweak sector of the standard model 
combines Coulomb and Higgs phases. Our current understanding of the phase structure of gauge 
theories owes much to the modern theory of phase transitions and critical phenomena, but has 
developed into a subject of extensive study. After reviewing some fundamental concepts of phase 
transitions and finite-temperature gauge theories, we discuss some recent work that broadly extends 
our knowledge of the mechanisms that determine the phase structure of gauge theories. A new class 
of models with a rich phase structure has been discovered, generalizing our understanding of the 
confinement-deconfinement transition in finite-temperature gauge theories. Models in this class have 
space-time topologies with one or more compact directions. On B? xS 1 , the addition of double-trace 
deformations or periodic adjoint fermions to a gauge theory can yield a confined phase in the region 
where the S 1 circumference L is small, so that the coupling constant is small, and semiclassical 
methods are applicable. In this region, Euclidean monopole solutions, which are constituents of 
finite-temperature instantons, play a crucial role in the calculation of a non-perturbative string 
tension. We review the techniques use to analyze this new class of models and the results obtained 
so far, as well as their application to finite-temperature phase structure, conformal phases of gauge 
theories and the large- N limit. 
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I. INTRODUCTION: GAUGE THEORIES 



Gauge theories are central to the modern understanding of the fundamental forces of 
nature, and form the foundation of the standard model of particle physics. Gauge theo- 
ries promote global symmetries, such as the U(l) invariance associated with conservation 
of electric charge, to local symmetries. In the case of quantum electrodynamics (QED), the 
the quantum field theory of electromagnetic interactions, local symmetry requires the intro- 
duction of a new vector field, associated with the photon. The standard model of particle 
physics is based on a local symmetry group of the form SU(2) x U(l) x SU(3). The stan- 
dard model contains the Weinberg-Salam-Glashow model of the electroweak interactions, 
with gauge group 577(2) x (7(1) , which unifies the electromagnetic interaction with the 
weak interaction. The number of gauge pariticles equals the rank of the gauge group, in 
this case four: the massless photon 7 and the massive vector boson W ± and Z°. The other 
part of the standard model is quantum chromodynamics (QCD), with gauge group SU(3), 
describes the strong interactions of quarks. There are eight gauge fields in QCD, describing 
the gluons that bind quarks inside of hadrons. It is widely believed that the standard model 
is not a complete description of particle physics, and there are many different proposals for 
what lies beyond the standard model. Gauge theories are integral to all such proposals, 
even for those where four-dimensional gauge theories are low-energy effective theories, as in 
string theory. 

One of the most fundamental questions we can ask about a given gauge theory is its 
phase structure. In the standard model, we observe three fundamentally different types 
of behavior: the familiar Coulomb behavior associated with the massless photon; the Higgs 
mechanism, responsible for the masses of the W ± and Z°; and the confinement of quarks and 
gluons by the gluon fields. These properties are characteristics of different phases: QCD 
is in a confined phase at zero temperature, while the electroweak sector of the standard 
model combines Coulomb and Higgs phases. The phase structure of gauge theories has been 
extensively studied. A large part of the interest in this subject centers on QCD at non-zero 
temperature and density, but there is now great interest in the phase structure of gauge 
theories associated with physics beyond the standard model. 

It is often convenient theoretically to view QCD in a simplified way, in a form without 
dynamical quarks; this is often referred to as pure S77(3) gauge theory. As a consequence 
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of dimensional transmutation, the dimensionless gauge coupling g 2 can be replaced by a 
parameter A with dimensions of energy. Key observables with units of mass such as hadron 
masses become pure numbers times A. In other words, the pure gauge theory is a theory 
with no adjustable dimensionless parameters, making it difficult to carry out analytical 
approximations. Including dynamical quarks with realistic masses does not make analytical 
work easier. However, there are many parameters that can be varied within the context of 
gauge theories, shedding light on important aspects of gauge theory behavior. Aside from 
their great experimental interest, finite-temperature gauge theories offer a dimensionless 
parameter T/A which can be used, for example, to change QCD from a confining phase 
at low temperatures to a plasma phase at high temperatures. This additional parameter 
tells us something about the physics of both phases. It is also natural to study the phase 
structure of gauge as other parameters, such as the number of colors N or the number of 
quark flavors Nf, are varied. 

Recent work has shown the existence of a new class of gauge theory models with important 
and desirable properties. All of these models have one or more compact directions, and the 
most developed case is the geometry R 3 x S 1 , which is the geometry of Euclidean gauge 
theories at finite temperature when the circumference L of S 1 is identified with the inverse 
temperature = 1/T. In addition to L, there are many parameters that can be used to 
examine the phase structure. Unlike conventional finite-temperature gauge theories, this 
new class can be put into a confined phase when L < A and g 2 (L) <C 1 [lj. Euclidean 
monopoles, the constituents of finite-temperature instantons, are essential to a semiclassical 
calculation of the string tension in this region. Moreover, this small-L phase is smoothly 
connected to the conventional, large- L confining phase |2]. There is a price to be paid for 
this. In the case of models with adjoint fermions, the use of periodic boundary conditions in 
the compact direction removes the spectral positivity of the transfer matrix in the compact 
direction, so L cannot be identified with /3 and must be regarded as a spatial direction. 
Models with double-trace deformations can be motivated as the heavy-quark limit of periodic 
adjoint fermions, and have a similar problem. The gains, however, are great: an analytic 
understanding of confinement in a class of four-dimensional models, the discovery of many 
new phases, and new approaches to conformality and to the large- A" limit. 

The analysis of this new class of models uses techniques and ideas from many areas 
of theoretical physics, and the results are broadly interesting as well. In section [TTj we 
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review the modern theory of phase transitions and critical phenomena, concentrating on 



key concepts and their appearance in a field-theoretic context. Section III introduces the 
basics of finite-temperature gauge theories, including important symmetries of gauge theories 



associated with confinement and chiral symmetry breaking. Sections IV through VI discuss 



recent developments. Section IV discussed the phase structure of gauge theories on R 3 x S 1 , 
an arena where it has proven possible to demonstrate confinement in certain models using 
semiclassical methods. SectionfVlfocuses on the general issue of gauge theory phase structure 



in the context of conformality and duality in gauge theories. In section VI, SU(N) gauge 
theories in the large- N limit are considered, and promising new approaches are compared 
with older formulations. A final section concludes and summarizes. 

The notations used throughout are as follows: All field theories are taken to be in Eu- 
clidean space unless noted otherwise. Lower-case Greek indices are used for space-time and 
the metric in Euclidean space is 

9^ = <r = V (!) 

Roman indices in the range j ■ ■ ■ n generally denote "spatial" directions on R 3 x S 1 , i.e., the 
three directions orthogonal to the compact direction. Roman indices in the range a - ■ ■ d 
generally label group generators, while capital letters are used to denote group representa- 
tions: F, Adj, S, A, R. S k is the fc-dimensional surface of a (k + l)-dimensional hypersphere, 
so S 1 is the unit circle. T k is the k- dimensional hypertorus, so T 1 is also S 1 . 



II. MODERN THEORY OF PHASE TRANSITIONS AND CRITICAL PHENOM- 
ENA 

A. Order parameters and symmetries 

The modern theory of phase transitions begins with Landau-Ginsburg theory j3] , and the 
key concepts of order parameters and free energies as functions of those order parameters. 
The simplest examples of order parameters come from ferromagnetic spin systems, where 
the magnetization rrij (x) is the order parameter, with j specifying the number of compo- 
nents. In order for the magnetization to be treated as a continuous variable, it is necessary 
to imagine that some smoothing procedure be applied to average the microscopic spins over 
many unit cells. The magnetization is naturally a three-component vector, but underly- 
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ing microscopic physics may cause it to be effectively one- or two-dimensional. Assuming 
macroscopic isotropy, the magnetization has a natural internal symmetry group 0(3), 0(2) 
or Z(2) depending on whether the order parameter is three, two or one dimensional. Thus 
the dimensionality of the order parameter need not be the same as the dimensionality of 
the system. Suppose for a moment we are interested in the free energy of a spin system 
in three spatial dimensions with a uniform magnetization which we take to be a scalar. 
This is the case where the magnetization has an easy axis; then <ft is the projection of the 
magnetization onto the easy axis. The key assumption of Landau- Ginsburg theory is that 
the free energy may be written as an integral over a local free energy density / ((f)) which is 
an analytic function of <fi as well as any underlying parameters, including the temperature. 
Thus the free energy is given by 

-4> + -< 

T 4 

where only terms even in (p are allowed by the Z(2) symmetry. If higher-order terms in the 
polynomial expansion are to be ignored, then we must have A > for stability. If r > 0, the 
free energy is minimized when = 0. As shown in figure [TJ if r < 0, the free energy will be 
minimized when 

* = ±7f P) 
giving the most common example of the idea of spontaneously broken symmetry. The free 
energy concept can be generalized to that of a function of a spatially-dependent magnetiza- 
tion <p (x). Including gradient terms in the expansion of the free energy, we obtain 



F[<j>] 



d 3 xf(0) = / d 3 x 



+ ... 



(2) 



d x 



(V0) 2 + 



+ 



A 



+ ... 



(4) 



where higher-order terms in V0 are suppressed. It is now possible to imagine calculating 
the free energy in an ensemble with an external field, via the functional integral 

d 3 x h (x) 4> (x) 



Z [h] = / [d(f>] exp 



+ 



(5) 



This is a three-dimensional Euclidean scalar field theory, the simplest example of the con- 
fluence of quantum field theory (at zero temperature) and classical statistical mechanics. 
Within Landau-Ginsburg theory, the point r = marks the location of a second-order phase 
transition, where the correlation length £ oc r" 1 / 2 becomes infinite. From a field-theoretic 
point of view, the inverse of the correlation length is a mass. From the modern point of 
view, a second-order phase transition is intrinsically associated with a massless excitation. 



f(0 

Figure 1: The Landau-Ginsburg free energy density f(4>) for the case r < where the symmetry 
(f) — > is spontaneously broken. 

B. Universality 

Landau-Ginsburg theory explains how non-analyticities arise in the vicinity of second- 
order critical point. For example, the critical point can be determined as a function of 
temperature: r (Tc) = locates the critical point. Analyticity of r in T gives r oc (T — T c ) 
near the critical point and hence (f) (T) oc \/T c — T in the low-temperature, symmetry-broken 
phase. The non-analytic behavior is captured in a set of critical indices, which characterize 
the critical behavior. For example, (T) oc (T c — T)' 3 so the critical index (3 is 1/2 in Landau- 
Ginsburg theory. Landau-Ginsburg theory leads to a simple universality: models with the 
same Landau-Ginsburg theory have the same phase diagram near the critical point, when 
parametrized in terms of the Landau-Ginsburg parameters, and the critical exponents at 
second-order phase transitions are idential. This universality was supported by the known 
exact equivalence between the Ising model and a lattice binary allow model, and led to 
the identification of the liquid-gas critical end point with the similar critical point in spin 
systems. However, both theoretical and experimental determinations of critical indices led 
to the conclusion that that critical indices, while falling into universal classes, were not given 
by the Landau-Ginsburg values [1]. Furthermore, it was clear that the spatial dimensionality 
d of the system must play a role. For example, classical spin systems with d = 1 with short- 
range interactions do not have phase transitions except at T = , in contradiction to the 
predictions of Landau-Ginsburg theory. It was clear that fluctuations, ignored in Landau- 
Ginsburg theory, were crucial. This led to the development of the modern renormalization 
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group [5]. 



C. Renormalization group 

The renormalization group takes many forms, but common to all approaches is the idea 
of a fixed point in the space of couplings. Consider an Ising model of spatial dimension d. 
It has two parameters: the lattice spacing a and a dimensionless coupling J which is the 
energy of a ferromagnetic bond divided by ksT. The Hamiltonian is given by a sum over 
nearest-neighbor interactions: 



where the spins Uj are associated with lattice sites j and take on the values ±1. On general 
grounds, the correlation length £ between spins must have the form £ = aF(J). At the 
critical point J c , the correlation length is infinite: F ( J c ) = oo. If an exact or approximate 
relation between £ and J is known, we can define a simple renormalization group transfor- 
mation as a change in the lattice spacing by a scale factor s > 1 such that a — > a' = sa 
and J — > J' (J) with the correlation length £ fixed: £' = £. In other words, we must have 



At the critical point J c , we must have J' = J, which is a fixed point of this renormalization 
group transformation. There are two equivalent points of view: From a physical point of 
view the correlation length is constant, but from a lattice point of view the dimensionless 
correlation length F(J), measured in units of lattice spacing, decreases away from J c . In 
the Ising model with d > 2, there are actually three fixed points. There are two fixed points 
with £ = 0, at J = and at J = oo, corresponding to T = and T = oo respectively. The 
second-order phase transition is associated with the non-trivial fixed point between them. 

It is generally necessary to introduce many parameters into the effective Hamiltonian 
describing the system at a length scale Xa, which can be denoted as the set {J a }- In Ising 
models, such additional couplings include a next-nearest-neighbor coupling and a four-spins 
around a square coupling. It is here that the renormalization group has its explanatory 
power. In the vicinity of a fixed point J* = {J*} , the renormalization group transform 




(6) 



rm 



sF (J') = F (J) or 



J> = F- X (s _1 F(J)) 



(7) 
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which changes the scale by a factor A can be linearized as 



(8) 



b 



The matrix R can be diagonalized and the corresponding linear combinations of spin inter- 
actions identified; the latter are usually referred to as operators and denoted O c . In the new 
basis for interactions the Hamiltonian can be written as "H = ^ c J c O c and a renormalization 
group transformation has the form J' c — J* = s Pc (J c — J*) . If p c > 0, then the renormaliza- 
tion group flow generated by J — > J' will take J c away from its fixed point value. In this 
case, the coupling and its corresponding operators are termed relevant. On the other hand, 
if p c < 0, the renormalization group flow will naturally take J c towards its fixed point value, 
at least within some basin of attraction where the linearization of the renormalization group 
is valid. Such couplings and operators are called irrelevant. There is also the case p c = 0, 
the marginal case, where higher-order effects in the renormalization group transformation 
must be considered. In the simplest case, a second-order phase transition is associated with 
the existence of one relevant coupling J r , accompanied by a potentially infinite set of ir- 
relevant couplings. This situation occurs, for example, in the Ising model in the absence 
of symmetry-breaking fields. Within the linear approximation, the hyper-plane in coupling 
space defined by J r = J*, is the critical surface where the correlation length is infinite. More 
generally, the set of renormalization group trajectories that flow into the fixed point forms a 
critical surface of codimension one in parameter space. The existence of irrelevant variables 
explains universality: many different models lie in the same universality class, and their 
critical behavior is determined not by the particular microscopic model, but by the physics 
near the fixed point. 

In the field-theoretic approach to critical phenomena, regularization of ultraviolet diver- 
gences typically requires the introduction of a parameter fi with dimensions of mass. Typical 
regularization schemes include a simple ultraviolet cutoff, Pauli-Villars regularization and 
dimensional regularization. The observed coupling constants are defined in terms of the 
arbitrary scale /i, and are therefore called running coupling constants. A change in \i does 
not change the underlying physics, but implies a corresponding change in the coupling con- 
stants, because a change in conventions does not change the physics of the system. In the 
case of a pure gauge theory, there is only one coupling g (//), and its running is described by 



dg_ 

dji 



(9) 
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where (3 (g) is the renormalization group beta function. In pure SU (N) gauge theories, 
the perturbative /3 function is negative, indicating that the renormalization group flow is 
towards g = at large /i. We say that g = is an ultraviolet fixed point, or that the theory 
is asymptotically free. There is an important connection here with the construction of the 
continuum limit of lattice gauge theories, where the inverse of lattice spacing a plays the role 
of an ultraviolet regulator. Physical quantities, such as a glueball mass M, cannot depend 
on the lattice spacing, but can be written parametrically as 



where g(a) is the coupling constant defined in the lattice action. The continuum limit, 
where M is to remain finite as a — > 0, can only be taken by simultaneously adjusting g 
such that F (g) — > 0. From the perspective of fixed lattice spacing, the continuum limit can 
only be taken at a second-order critical point, and physical masses become small in lattice 
units. Although the most common gauge theories have a simple fixed point structure in four 
dimensions, non-trivial fixed points play an important role in some scenarios for physics 
beyond the standard model, as will be discussed in section |V} 

The close connection between quantum field theories at zero temperature in d Euclidean 
dimensions with classical statistical mechanics models was crucial in developing the field- 
theoretic approach to critical phenomena. The simplest example of this connection is be- 
tween the simple model of ferromagnetism described above and the <ft A field theory. There is 
a simple rule of thumb which is helpful in matching up field-theoretic properties with renor- 
malization group properties: a super-renormalizable interaction is relevant, a renormalizable 
interaction is marginal, and a non-renormalizable interaction is irrelevant. For example, in 
the of the the A0 4 field theory in d = 4 — e dimensions |5HE], the 4 interaction is marginal 
in four dimensions, with an infrared fixed point at A = 0. For e > 0, a non-trivial infrared 
fixed point A* > is present, and the (ft 2 term is the only relevant operator near the critical 
point. From this point of view, Landau-Ginsburg theory fails to give the correct critical 
indices for d < 4 because it is essentially an expansion around A = rather than A*. Above 
four dimensions, the IR fixed point is at A = 0, and the critical indices are correctly pre- 
dicted by Landau-Ginsuburg theory. The field theoretic approach allows a precise statement 
of universality: starting from the most general renormalizable field theory using the order 
parameter as the field, the second-order phase transitions and their properties follow from 




(10) 
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the fixed points of the theory. Thus, the universality class of a second-order phase transition 
depends on the dimensionality of space and the symmetry group of the order parameter. 

Although the renormalization group has led to tremendous progress in our understanding 
of second-order phase transitions, the vast majoritiy of phase transitions observed in nature 
are first-order. Such transitions are generally accompanied by a discontinuity in one or 
more order parameter, and there is a latent heat associated with the transition. Within the 
framework of the effective potential, first-order transitions occur when there are two or more 
disconnected global minima of the free energy. A simple model that illustrates the general 
concept has a Landau- Ginsburg free energy density of the form 

/(0) = ^ 2 + ^ 4 + |0 6 (11) 

Unlike the simpler 4 theory, with g > to ensure stability we need not have A > 0. For 
A > 0, there is a second-order phase transition at r = 0. However, if r < 0, there is a 
first-order phase transition at A = — A^Jrg/3 where two non-trivial minima have the same 
free energy as the local minimum at = 0. 

This model contains an important lesson: Although the behavior of systems at second- 
order transitions is determined by the local structure of the the free energy density /(</>), 
the overall phase structure is determined by the location of the global minimum of the free 
energy. This is usually more difficult to obtain in systems with several order parameters and 
complicated interactions. 



D. The effective potential 



The close correspondence between classical statistical mechanics and quantum field theory 
extends beyond critical phenomena. In thermodynamics and statistical mechanics, different 
independent variables, free energies and statistical ensembles are often used. These ideas, 
so useful in statistical physics, have direct analogs in quantum field theory. We begin with 
the case of quantum field theories at zero temperature, and later extend our results to finite 
temperature. Consider the following generalization of the partition function 



Z [J] = / [dtp] exp 



-S 



d d+l x J (x) (j) (x) 



(12) 
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This is a generating functional for all n-point correlation functions. 

fjif = 1 + E^[/ d d+1 x 1 ..d d+1 x n J(x 1 )..J(x n )(<f>(x 1 )..<P(x n )). (13) 

The logarithm of Z [J], —W [J] = \ogZ [J] is the generating functional of connected corre- 
lation functions 



00 1 r 

-W[J\ = -W[0]+^2~ / d d+1 Xl ..d d+1 x n J (xi) ..J (x n ) (0 (an) ..0 (x n )) c (14) 

n=l 

where the connected correlation functions are, for example, 

(0 M (x 2 )) c = (0 (xx) (x 2 )> - (0 (xx)) (0 (x 2 )) (15) 

While [0] is typically not physically meaningful, the difference W [ J] — W [0] can be 
interpreted in the case of constant J as VTS(J), where VT is the volume of Euclidean 
space-time, and £ (J) is the vacuum energy density in the presence of the constant source 
J. In this formalism, J (x) plays a role analogous to that of a spatially dependent external 
magnetic field in a ferromagnet, and (0(x)) is the analog of the local magnetization at x. 
Spontaneous symmetry breaking is typically determined from the solution of the relation 

SW 

<*<*» = -iTW™ (16) 

This is generally not a closed equation for (0(x)), and approximation solutions from, e.g., 
perturbation theory must be used. The functional W [J] is the field-theoretic analog of the 
free energy of a spin system in an external magnetic field. It will be convenient in what 
follows to work with as the independent variable rather than J. The effective action T [0 C ] 
is a functional of the classical field C (x) defined by 

5W 

*<*> = -«M (17) 

for arbitrary J (x). Given that relation between J(x) and (f> c (x), we define T via a functional 
Legendre transform: 

T [0 C ] = W [ J] + j d d+1 x J (x) C (x) (18) 

which satisfies 

m - J{x) - (19) 
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This construction is precisely analogous to exchanging the role of independent variable be- 
tween the external field and the magnetization in a ferro magnet. For constant fields, W[J] 
is the free energy of the system with constant J, and /Gamma[(f)} is the free energy of the 
system with the conjugate variable <p c held constant. For the case of a single scalar field, 
the effective potential is given at one loop as 

T [0 C ] = S [0 C ] + \tt [log (-V 2 + V" (&))] (20) 

r also has a graphical interpretation as the generator of one-particle- irreducible (1PI) Feyn- 
man graphs. 

In many applications, it is convenient to consider the effective potential V e ff(<f> c ). It can 
be obtained from a derivative expansion of T[(p c ] for slowly- varying fields: 



/ 



d d+1 x 



V eff (<j> c ) + ^Z(<j> c ) (d^ c ) 2 + ... 



(21) 



where quartic-derivative terms and higher are indicated by ellipses. The effective potential 
should be understood at zero temperature as the vacuum energy density of the system given 
that (4>(x)) = 4> c . We have for constant fields 

J (22) 



dV eff 



d<p c 

and V e //(0 C ) is related to £(J) by a Legendre transform 

V eff ((j> c ) = £{ J) + J<f> c . (23) 

It is important to note that this construction is ambiguous when there are multiple solutions 
for C . In general, V e //(0 C ) is given by the convex hull of the differential construction, a 
generalization of Maxwell's equal construction. 

V e ff is an extension of the classical potential V including quantum effects. From the 
point of view of perturbation theory, V e ff = V + 0(h), and the 0(h) correction term has .... 

v eff (&) = V (0 C ) + \ J 0- 4 log [k 2 + V" (0 C )) (24) 

Note that this expression requires renormalization. It is easy to show that this expression 
may also be written as 

V eff (0 C ) = V (0 C ) + hj Vk 2 + V" (0 C ) (25) 

which clearly exposes the one-loop correction as a sum over the zero-point energies of all the 
field modes. The background field method offers a very powerful technique for evaluating 
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E. Field theory at finite temperature 



In quantum statistical mechanics, the partition function is given by Z = Trfexp (—/3H)] 
where /3 is the inverse of the temperature T. The propagation in Euclidean time is over a 
finite extent, given by /3, and the trace can be implemented by requiring periodic boundary 
conditions on the fields. The partition function is just the generating function of the field 
theory 

Z = J [d(p}e- s (26) 

but the functional integral is now over all field configurations satisfying <f> (x, 0) = <f>(x,/3). 
Thus the effects of non-zero temperature are included by changing the geometry of the 
system from R d+l to R d x S 1 with periodic boundary conditions for bosons; fermions require 
antiperiodic boundary conditions. 

The generalization of the one-loop effective potential to finite temperature is illuminating. 
Due to the periodic boundary conditions in the Euclidean time direction, the integral over 
is replaced by a sum over Matsubara frequencies 

K//(0 c ) = m) + ~£/ ^ h ^[{ 2 f) 2 + p + v ''^) W 

which can be converted into the form 

d 3 k 



/d k I P» / 
-V^ + ^"(0 c )+log(l-e 



-Py/k 2 +V"(<t> c ) 



(28) 



(2tt)< 

The additional term is instantly recognizable as —p, where p is the pressure of a relativistic 
Bose gas with mass VV"(0 c ). This is of course consistent with the identification of V e ff (0 C ) 
with the free energy of the system, when evaluated at its global minimum. 



F. Symmetry restoration at high temperatures 

Generally speaking, broken symmetries are associated with low temperatures, as seen, 
for example, in spin systems such as the Ising model. This can be explicitly seen in scalar 
field theories at finite temperature using the effective potential. The finite temperature part 
of the one-loop effective potential has a natural expansion for high temperatures. A more 



general form for such expressions will be given in section III C so for now we simply note 
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that that for a single real scalar field the first two terms in a high-temperature expansion 
are given by P, [TU] 

Vi ^- w -w + h T2y " (4 " ) - (29) 

The first term is the usual blackbody term with no effect on the value of C . However, the 
sub-leading T 2 term acts to restore symmetry. If the potential V(0) has the form 

^(0) = -^V + ^A0 4 , (30) 

then /i 2 > will lead to spontaneous symmetry breaking at zero temperature. However Vyr 
supplies a temperature-dependent mass term such that the potential for becomes 

v (0) + v 1T (0) - + 1 (At 2 - ^ 2 + I A 4 (si) 

and the symmetry — > —0 will be restored at temperatures above T c ~ ^/24/A/x. A 
notable exception to the general rule of symmetry restoration at high temperature is the 
deconfinement transition in pure gauge theories, where the broken phase lies above the 



unbroken phase in temperature. This will be explain in detail in Section HI B 



The temperature dependence of the effective potential is very similar to the general form 
of the free energy as a function of order parameter postulated by Landau theory. A slightly 
different construction, often referred to as dimensional reduction, also leads to a free energy 
of Landau's form. In dimensional reduction, the field modes with non-zero Matsubara 
frequencies are integrated out, so that only modes which are constant in Euclidean time 
remain. This effectively reduces a (d + l)-dimensional theory to a d- dimensional theory 
with temperature-dependent parameters in the action. 



G. Spatially modulated phases 

A further complication in the determination of the phase diagram is the possible existence 
of spatially modulated phases. In the most familiar cases, symmetry breaking occurs when 
some order parameter (x) has a constant, non-zero expectation value independent of x. 
However, there are phases in nature where order parameters are not translationally and/or 
rotationally invariant. The most familiar example is the formation of crystalline solid phases 
in ordinary materials, where the order parameter is the density. In a liquid or vapor phase, 
the density is uniform, but in a crystalline phase it is periodic. Magnetic systems can also 
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show spatially modulated phases, as exhibited by simple models such as the anisotropic 
next-nearest-neighbor Ising (ANNNI) model [11J and the chiral potts moder [T2"j [TB"] . In 
field theories, spatially modulated phases are generally associated with systems at finite 
density, and the sign problem is intimately involved in this behavior [bfl [15]. In QCD, 
some color superconducting phases are crystalline in nature [IE]; these are analogs of the 
Larkin-Ovchinnikov-Fulde-Ferrell (LOFF) phases of ordinary superconductors [17, 18j. 

There is a simple approach to spatially modulated phases due to Lifshitz [3j. As a simple 
model, consider a Landau-Ginsburg model with a free energy density of the form 

/(0) = ^W) 2 + ^(VV) 2 + ^ 2 + ^ (32) 

where b and A must be positive but a and r can be of either sign. Taking the case a < and 
r > 0, it is easy to see that the homogeneous equilibrium state will be unstable to spatially- 
varying perturbations of wavenumber k for any values of k satisfying bk 4 + ak 2 + r < 0. This 
model has three phases: an unbroken phase, a broken phase, and a spatially modulated 
phase. The point a = 0, r = is the Lifshitz point, where the three phases coexist. 

H. Non-equilibrium behavior: nucleation, spinodal decomposition and relaxation 

The perturbative effective potential can yield information about non-equilibrium behav- 
ior. Recall that the perturbative effective potential does not satisfy the convexity properties 
required of V e ff by its Legendre transform construction. However, the non-convex region of 
V e ff contains information about non-equilibrium behavior. Metastable states can be under- 
stood as local minima of the effective potential. Consider a field theory with a double well 
potential of the form V ((f)) = X ((f) 2 — v 2 ) 2 . The presence of a small additional linear coupling 
— J(f> will bias the system towards one of the two minima at ±v . This leads to a discontinuity 
in the behavior of C as a function of J. However, the perturbative effective potential shows 
the persistence of metastable states across J = 0. The intuitive picture is simple in the 
thin wall approximation, which has roots in classical nucleation theory. Consider an initial, 
homogeneous field configuration where = (f>i OC ai, the value of the local minimum of V e ff. 
Coherent fluctuations in space will produce small droplets, also called bubbles, of (f) g i bai, the 
global minimum. The difference between the free energy densities Af = f((f>i OC ai) — f (<f> global) 
drives the expansion of such bubbles. However, there is also a surface tension a residing in 
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the interface between the two phases, which opposes the growth of bubbles. In the thin-wall 
approximation, the free energy of a three-dimensional droplet of radius R is approximately 

Air 

F (R) = -Af ■ —R 3 + a ■ 4rf (33) 

where a is the interfacial surface tension at equilibrium. The surface tension may be cal- 
culated from the one- dimensional kink solution that interpolates between the two phases 
A droplet which is smaller than the critical size R c will shrink due to surface tension ef- 
fects, while a droplet with R > R c will expand to convert all space to the the equilibrium 
phase. The decay rate of metastable states can be calculated within the functional integral 
formalism using instanton techniques [X9H2T] . 

r = A e~ Sb (34) 

where Sj, is the action of the so-called bounce solution. Within the thin-wall approxima- 
tion, it is equal to F(R C ). The pre- factor A is given in terms of a functional determinant 
representing fluctuations around the bounce solution. 

Within the saddle-point approximation considered here, there is a clear point at which 
the metastable phase ceases to exist, i.e., where the metastable states cease to exist as 
local minima. The subspace of parameter space where this occurs is referred to as the spin- 
odal, which may be a point, line, et cetera. The spinodal can be located by the condition 
V eff (4>s P inodai) = V" ff ((f) spinodal) = 0, and more generally by the vanishing of a mass asso- 
ciated with a metastable state. In physical systems, fluctuations become important near 
the spinodal point, and the distinction between metastable versus unstable states becomes 
blurred. 

Unstable states are of interest in their own right, and something about their behavior 
can be learned using the effective action. Consider a 4 field theory in d spatial dimensions 
at finite temperature. We apply dimensional reduction so we have a three-dimensional field 
theory with temperature-dependent coefficients in the potential. Suppose the system begins 
in a high-temperature equilibrium state with unbroken symmetry so <fi = 0, but has a stable 
broken symmetry phase at low temperatures. A rapid lowering of the temperature below 
the critical point, known as a rapid quench, will place the system in a state of unstable 
equilibrium, with V"tt (0 = 0) < 0. As in the case of a metastable state, such an unstable 
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state will lead to a free energy 

f(<j ) = 0) = V eff (<j ) = 0) = V( ( j ) = 0)+ 1 - I -0^1og(P + V''(0 = O)) (35) 

— * 

Because V"jf (0 = 0) < , the modes near k = in the functional determinant give rise to 
an imaginary component of the free energy This imaginary part can be associated with the 
decay rate of the unstable phase in a manner similar to the decay of metastable states [22J . 
In the early stages of equilibration, the unstable modes in the region k 2 + V" (<p = 0) < 
lead to exponentially growth of long-wavelength oscillations in the order parameter. This 
process is known as spinodal decomposition. When the order parameter is not conserved, 
e.g. in the case of magnetization, the fastest growth occurs around k = 0; however, in the 
case of a conserved order parameters, such as a charge density, the fastest growth occurs at 



a non-zero value of 



m 

In addition to nucleation and spinodal decomposition, there is a third equilibration pro- 
cess, one which is always present in non-equilibrium situations. Relaxation is the process 
whereby a system returns to its equilibrium value after a small perturbation away from 
equilibrium; the term may also be applied to individual field modes. An introductory dis- 
cussion of relaxation requires a treatment of damping, i.e. frictional or dissipative terms in 
the effective action. All three processes have been extensively discussed in the context of 
the early universe [231 [21], particularly with respect to inflation. 



III. PHASES OF FINITE TEMPERATURE GAUGE THEORIES 

Gauge theories with T ^ ("finite temperature") have a rich phase structure which has 
been extensively explored using a combination of analytic methods and lattice simulations. 
Non-Abelian gauge theories have global symmetries and associated order parameters which 
are analogous to magnetization in spin systems, and much of the modern formalism of 
critical phenomena is directly applicable. The symmetries and order parameters associated 
with quark confinement and chiral symmetry breaking are of particular interest as principal 
determinants of gauge theory phase structure. 
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A. Confinement 



The most striking feature of QCD is confinement: the force between widely separated 
quark- ant iquark pairs is a constant a, known as the string tension. This constant force 
implies a potential energy between a quark- ant iquark pair that grows as or for large distances 
r. The numerical value of a, determined from phenomenology, is approximately 0.18 GeV 2 ~ 
0.9 GeV/ fm. The string tension can be extracted from the asymptotic behavior of the 
quark- ant iquark potential. In turn, the static potential between two heavy quarks can be 
determined in a gauge-invariant way from the expectation value of the Wilson loop. This is 
a non-local operator associated with a closed curve C in space-time parametrized as x ^ 
where r can be taken to be the unit interval and x M (1) = x M (0). The Wilson loop is defined 
as 



where V indicates path-ordering of the gauge fields (x) along the path C |25| . Taking 
W [C] to be an element of the abstract gauge group G, the basic observables are TtrW [C] , 
where the trace is taken over an irreducible representation R of G. Tr^W [C] has a physical 
interpretation as the non-Abelian phase factor associated with a heavy particle in the rep- 
resentation R moving adiabatically around the closed loop C. Typically, we are interested 
in rectangular Wilson loops of width sides L and T. The loop can be associated with a 
process in which a particle-antiparticle pair are created at one time, move to a separation 
L, propagate forward in time for an interval T, and then annihilate. For T ^> L, we have 



where Vr (L) is the heavy quark-antiquark potential for particles in the representation R. 
For large distances, the potential can grow no faster than linearly with L, and that linear 
growth defines the string tension gr for the representation R: 



as L — » oo. The O (1) correction term represents a perimeter-law contribution, which must 
be renormalized; the effect is closely related to the mass renormalization of a heavy particle 
interacting with the gauge field. 

Beyond the perimeter law contribution, there are corrections which are 0(1/ L) and 
higher that are of physical interest. Using an effective string description, Luscher showed 




(36) 



(Tr R W[C])~exp[-V R (L)T\ 



(37) 



V (L) ->■ o R L + O (1) 



(38) 
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that the O (1/L) has a universal coefficient that depends only on the dimensionality of space- 
time |2S]; this prediction was confirmed in lattice simulations |27j . Subsequent analytical 
work has explored the effects of higher-order terms |2"%H3"0] . There have also been results on 
the broadening of the flux tube as a function of L [3T| 132] . For a lattice-oriented review of 
string effects in gauge theories, see the review by Pepe |33| . 

The statement that quarks are confined is the statement that the string tension op of the 
fundamental representation of SU (3) gauge theory is non-zero. However, each representation 
has its own string tension. In lattice simulations of S77(3), the fundamental string tension is 
the smallest, and representations of larger dimension generally have a larger string tension. A 
characteristic feature of gauge theories in four dimensions is dimensional transmutation: the 
classical action is scale-invariant, but the introduction of a mass scale // is necessary in order 
to define the running coupling constant <?(//). It follows directly from the renormalization 
group that any physical mass in four-dimensional pure gauge theories must be proportional 
to a renormalization group invariant mass, usually written simply as A. This implies that 
a R — CrA 2 , where cr is a pure number. Put another way, the ratio Ur/uf is a pure 
number characterizing the representation R. In (1 + 1) -dimensional gauge theories, the string 
tension or is proportional to the quadratic Casimir invariant, leading to a behavior known 
as Casimir scaling. It is generally believed that a satisfactory description of confinement in 
four dimensions would include an understanding of the string-tension scaling law. 

There are two subtleties associated with the string tension. The first subtlety is that 
string tensions are generally determined in so-called pure gauge theories, where only the 
gauge fields are dynamical; quarks and other particles exist only as static classical charges. 
This is done because of string breaking , also known as charge screening. Consider a theory 
with dynamical particles of mass m in the fundamental representation F of the gauge group. 
Then energy obtained by separating a pair of static sources in the fundamental representation 
is opL\ when that energy becomes on the order of 2m, it will become energetically favorable 
to produce a pair of dynamical particles at a cost in energy of 2m, and no string tension 
will be seen for L > 2m ja p- 

The second subtlety is related: in SU (N) gauge theories, a state with N + 1 fermions in 
the fundamental representation, e.g. quarks in 577(3), will form a baryon, a neutral state of 
N fermions, plus one additional fermion. Thus the string tension between N+l fermions and 
N + 1 antifermions observed at large distances should be the same string tension observed 
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between one fermion and one antifermion. Thus for purposes of confinement, the relevant 
"charge" is an integer k which can be taken between 1 and N. Furthermore, the string 
tension is the same for states of iV-ality k and N — k, so that there are [yl independent 
string tensions. This behavior is associated with the Z(N) symmetry underlying confinement 
in SU(N), discussed below. 

It is generally believed that there is some simple regular behavior described by string 
tension scaling laws [31]. For example, in (1 + l)-dimensions, Casimir scaling gives 

a k _ k(N-k) 

a t ~ N-l [ } 

This behavior is also found in the strong-coupling limit of lattice gauge theories in all 
dimensions, and is consistent with the behavior seen in lattice simulations. An alternative 
to Casimir scaling is sine-law scaling 

"± = ^3) (40) 
<x, sin(f) 

a behavior predicted in certain supersymmetric gauge theories [35], in MQCD [36]. This 
behavior is also roughly consistent with lattice results. Note that both formulas give cr/v = 
and are invariant under k -h- iV — k. In the limit N — > oo at fixed k, both formulas give 
<7fc — » ko\. Corrections to the large- N limit start at 1/N for Casimir scaling and 1/N 2 for 
sine-law scaling, and it has been argued that the large- N limit requires 1/N 2 corrections 
|37| . However, lattice simulations of SU(N) gauge theories in (2 + 1) dimensions indicate 
that 1/N corrections are strongly indicated |38j . and a careful analysis shows that Casimir 
scaling is compatible with the large- N limit 



B. Pure gauge theories at finite temperature 

If one or more directions in space-time are compact, the string tension may also be 
measured using the Polyakov loop P, also known as the Wilson line. The Polyakov loop 
is essentially a Wilson loop that uses a compact direction in space-time to close the curve 
using a topologically non-trivial path in space time, as shown in figure [2j The typical use for 
the Polyakov loop is for gauge theories at finite temperature, where space-time is R 3 x S 1 . 
The partition function being given by Z = Tr [e - ^] , the circumference of S* 1 is given by 
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Figure 2: The Polyakov loop is associated with the worldline of a heavy particle. 



the inverse temperature = 1/T. In this case, we write 



P (x) — V exp i / dx^A^ (x) 

L Jo 



(41) 



and string tensions may be determined from a two-point function 



(Tr R P (x)Tr R P^ (y)) 



~ e 



(42) 



a behavior that assumes that the one-point function (TtrP (x)) = 0. This behavior is 
shown in figure [3j The Polyakov loop one-point function (TtrP (if)) can be interpreted as a 
Boltzmann factor exp (—/3Fr), where Fr is the free energy required to add a static particle 
in the representation R to the system. Of course, (TtrP (x)) = implies that Fr = oo , 
which is thus a fundamental criterion determining whether particles in the representation R 
are confined. Note that the introduction of a compact direction breaks the four-dimensional 
symmetry of the theory, and the string tension measured by Polyakov loops is not the same 
as the string tension measured by Wilson loops lying in non-compact planes. In the case of 
finite temperature, it is natural to use the terminology electric and magnetic string tension, 
respectively. In the limit where the compactification radius becomes large, e.g. (5 — > oo, the 
two string tensions must coincide. 

1. Center symmetry 

One of the most important concepts in our understanding of confinement is the role of 
center symmetry. The center of a Lie group is the set of all elements that commute with 
every other element. For SU(N), this is Z(N). Although the Z(N) symmetry of SU(N) 
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Figure 3: The Polyakov loop two-point function determines the electric string tension. 

gauge theories can be understood from the continuum theory, it is easier to understand from 
a lattice point of view. A lattice gauge theory associates link variable (x) with each lattice 
site x and direction //. The link variable is considered to be the path-ordered exponential 
of the gauge field from x to x + fx: (x) = exp [iaA^ (x)] . Consider a center symmetry 
transformation on all the links in a given direction on a fixed hyperplane perpendicular to 
the direction. The standard example from SU(N) gauge theories at finite temperature is 
U4 (x,t) — > zUi (x,t) for all x and fixed t, with z G Z(N). Because lattice actions such as 
the Wilson action consist of sums of small Wilson loops, they are invariant under this global 
symmetry. However, the Polyakov loop transforms as P (x) — > zP(x), and more generally 

Tr R P (x) z kR Tr R P (x) (43) 

where k R is an integer in the set {0, 1, N — 1} and is known as the iV-ality of the rep- 
resentation R. If kft 7^ 0, then unbroken global Z(N) symmetry implies (Tr R P (x)) = 0. 
Thus global Z(N) symmetry defines the confining phase of a gauge theory. For pure gauge 
theories at non-zero temperature, the deconfinement phase transition is associated with the 
loss of Z(N) symmetry at the critical point T^. Below that point (TrpP (x)) = but above 
T d , (Tr R P{x))^0. 

Notice that the case of zero iV-ality representations is special within this framework: 
there is no requirement from Z(N) symmetry that these representations are confined. This 
includes the adjoint representation, the representation of the gauge particles. However, 
lattice simulation indicate that (Tr R P(x)) is very small for these representations in the 
confined phase. Although screening by gauge particles must dominate at large distances, 
these zero iV-ality representations have well-defined string tensions at intermediate distances 
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scales, e.g., on the order of a few fermi for SU(3), behaving in a manner very similar to 
representations with non-zero iV-ality [40, 41 J . 



2. The nature of the confined phase and partially confined phases 

The requirements for confinement are simple: Z(N) symmetry and a mass gap. However, 
gauge theories with SU(N) symmetry are not simply Z(N) systems. In models where the 
order parameter is an element of Z(N), the order parameter must average to zero in the 
disordered phase, and the equilibrium disordered state must consist of many domains, each 
with its own Z(N) orientation. This need not be the case with SU(N) systems. There is a 
unique set of Polyakov loop eigenvalues that are Z(N) symmetric, and can form the basis 
for a semiclassical understanding of confinement. 

Although TrpP is the order parameter for confinement of particles in the fundamental 
representation, it does not by itself characterize the confined phase, and the most general 
form of the free energy does not depend on P solely through TrpP and TrpP^ [12] ■ The 
high-temperature effective potential V e ff (P) illustrates this point. It is not a function solely 
of TrpP and its conjugate, and cannot be written as an infinite series in TrpP and TrpP^ . 
A simple example will illustrate this point. Consider two diagonal matrices lying in S77(4), 
defined by 

/ e i7r/4 \ 



Pi 



3 i37r/4 



i5n/4 



(44) 



and 



3 i7r/2 



\ 



-itt/2 



V 



-in/2 



(45) 



Both Pi and P2 have zero trace in the fundamental representation, yet the traces of their 
square are differenct: TrpP 2 = and TrpP 2 = —4, establishing that TrpP 2 cannot be a 
function solely of TrpP. 
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At first sight, this may seem to contradict two standard results: a) the characters form 
a complete, in fact orthogonal, basis on class functions; b) all characters may be obtained 
from the fundamental representation by repeated multiplication and the application of 

Xa (P)x b (P) = ]Tn(a,&;c)x c (P) (46) 

c 

where all n's are non-negative integers. Taken together, these results might suggest that all 
group characters are polynomials in Tr F U and its complex conjugate. Consider, however, 
the product representation N®N. It is reducible into N(N + l)/2 © N(N - 1) /2, which are 
symmetric and antisymmetric representations, respectively. Their characters are respectively 

Xs(P) = \[{Tr F P) 2 + Tr F {P 2 )] 

Xa(P) = \[{Tr F P) 2 -Tr F {P 2 )]. (47) 

Note that the sum xs + Xa is a polynomial in Tr F P, but xs an d Xa are in general not. For 
example, the restriction on the eigenvalues of P imposed by det (P) = 1 allows us to prove 
that for 51/(3) 

l - [{Tr F Pf - Tr F (P 2 )] = Tr F (pt) (48) 

in accord with the SU(3) result 3 © 3 = 6 © 3. In this case, it is true that Tr F (P 2 ) can 
be written as a polynomial in Tr F P and Tr F P* . However, for S77(4), unitarity of P gives 
instead 

I [{Tr F Pf - Tr F (P 2 )] = \ [{Tr F Pf - Tr F (P 2 )]* (49) 

which shows that the 6 representation of SU(A) is real. SU(N) group characters can be 
written as polynomials in Tr F (P) and its complex conjugate only for 577(2) and 577(3). 

An alternative statement is that the Polyakov loop in the fundamental representation, 
Tr F U, is not sufficient to determine the eigenvalues of P, beginning with the case of SU(A). 
Let us label the eigenvalues for the Polyakov loop as z i: with % — 1 to N. For SU(3), the 
characteristic polynomial for the Polyakov loop is 

3 

Y\(Z - Zi) = z 3 - (zi + z 2 + z 3 ) z 2 + (z ± z 2 + z 2 z 3 + z 3 z ± ) z - z 1 z 2 z 3 . (50) 
1=1 

Since the determinant of a special unitary matrix is 1, we have -21^2-23 = 1, and the charac- 
teristic polynomial is 

z 3 - z 2 Tr F P + zTr F P + - 1 (51) 
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so for SU(3), knowledge of TtfP determines all eigenvalues, and the free energy can be writ- 
ten as a function of TrpP alone. For SU(4), similar considerations allow the characteristic 
polynomial to be written as 



and knowledge of TrpP must be supplemented by the value of Trp (P 2 ). As N increases, 
more information must be supplied to reconstruct the eigenvalues. 

The characteristic equation for Polyakov loops in SU(N) may be written as 

N 



where Xj is the character for the representation formed from the antisymmetrized product of 
the fundamental j times; the corresponding Young tableaux is j vertical boxes. Of course, 
Xj (P) has iV-ality j, and Xo{P) — Xn (P) = 1- At the classical level, a Polyakov loop 
that represents the confining phase must have Xj (P) = for 1 < j < N — 1, leading 
to a characteristic equation 1 + (— A)^ = 0. This in turn implies the important result 
that any model in which confinement is associated with a particular field configuration (up 
to gauge transformations) must have Polyakov loop eigenvalues evenly spaced on the unit 
circle in the confined phase. This is precisely what must happen in the large-iV limit, and 
is known to occur in analytically tractable lattice models of the large- N limit |43j . This 
behavior is not possible in simple Z(N) models where symmetric phases always represent 
some average over different Z(N) domains. However, it should be noted that there are Z(N) 
spin systems with features similar to models with variables in SU(N). The best-known is 
the Blume- Emory- Griffiths model, which was originally proposed as a simple model of 3 He- 
A He mixtures |44| . Essentially, an Ising spin model is extended such that each spin a can 
have the value as well as ±1. The behavior of such extended Z(N) lattice models is very 
similar to treatments of SU(N) lattice models with fundamental and adjoint interactions, 
where the adjoint interaction can be adjusted to give a Z(N) limit |43| 136] ; see also |47j . 

In the confining phase of a gauge theory, Z(N) symmetry is unbroken, and all repre- 
sentations with non-zero A r -ality are confined. In the deconfined phase, Z(N) symmetry is 
completely lost, and particles are no longer confined, independent of their representation. 
For N > 4, additional phases are possible where Z(N) is broken down to a non-trivial sub- 
group [2, 42\ HHl H9]. In the case of Z(4), there can be breaking of center symmetry down to 




(52) 




(53) 
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Z(2). In this partially confined phase, states consisting of two fundamental representation 
fermions are not confined, but single fermions are. In the case of SU(N), Z(N) can break 
to Z(k), where k is any divisor of N. States with 1 to k — 1 fermions are confined, but states 
with k fermions are not. It is often convenient to include the confined and deconfined phases 



as k = N and k — 1, respectively. As we will see in section [TVj such partially confined phases 
been found in gauge theories on R 3 x S 1 , using both lattice simulations and perturbation 
theory [2J. It should be noted that not all gauge groups have non-trivial centers. The gauge 
group G(2) provides an interesting example of a gauge theory without a center |5Uj that has 
received significant attention [ST]. 



3. Origin of the deconfinement transition 



As we have seen in section [TTJ the critical properties of a four-dimensional field theory at 
finite temperature can often best be understood in terms of an effective three-dimensional 
theory. From this point of view, the Polyakov loop can be considered to be a scalar order 
parameter with Z(N) invariance in some effective three-dimensional theory. A complete 
treatment of the effective potential, including non-perturbative effects, would presumably 
yield the phase structure, including both the confined and deconfined phases. 

Perturbation theory is a reliable indicator of broken center symmetry and thus deconfine- 
ment at high temperature, because the running coupling constant g(T) is small if T 3> A. 
The one-loop effective potential for a pure gauge theory in the background of a static 
Polyakov loop P can be easily evaluated in a gauge where A4 is time-independent and 



diagonal j52j [53]. It is easy to see that V})t (P) is given by 



v eff (P)=2TTr Adj J2 



d 3 k 



log 



(27mT - A A f + k 2 



(54) 



nez" (2*0 

where the factor of 2 represents the two helicity states of each mode. Note that there is 
no classical contribution. Discarding the zero-point energy term, we obtain the one-loop 
finite-temperature effective potential for gauge bosons 

d 3 k 



V e 1 } f (P) = 2TTr 



Adj 



(27T) 



log 



1 — P exp 



/T 



(55) 



which is the free energy density of the gauge bosons in the background of the Polyakov loop 
P. 
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The logarithm in this expression for V^jj (P) can be expanded, leading to an interpretation 
of Vgjcy (P) as a sum of contributions from gluon worldlines wrapping around the compact 
direction an arbitrary number of times. Explicitly, we have the expression 



2 ^ 1 



v eif(p)--- 2 }2-J^ pn - ( 56 ) 

n=l 

From this form, it is easy to see that V^jf (P) is minimized when all the moments Tr^P n 
are maximized. This occurs when P G Z(N), which gives Tr^P™ = N 2 — 1. This indicates 
that the one-loop gluon effective potential favors the deconfined phase. The pressure p is 
the negative of the free energy density at the minimum, 

P (T) = 2(iV 2 -l)^, (57) 

which is exactly p for a blackbody with 2 (iV 2 — 1) degrees of freedom. 

In the gauge where A4 is diagonal and time-independent, we can parametrize A4 in the 
fundamental representation of SU (N) as a diagonal, traceless N x N matrix 

(A 4 ) jk = T6 3 5 jk (58) 

so that 

Pjk = e id i5 jk (59) 

with YljLi @j = 0- Using the decomposition F (g) F = 1 © Adj and the corresponding 
decomposition 

Tr F PTrP+ = l + Tr Adj P, (60) 

one realizes that the iV 2 eigenvalues of the product representation F (g) F have the form 
exp [iA9j k ] , where we define A9j k = 9j — k . The effective potential is given by 

N 00 

V eff (P) - ~~ 2 E (! - N 6 ^ E ^ 6X P [ <nA ^l • ( 61 ) 
j„fc=l n=l 

The infinite sum over n may be carried out explicitly in terms of the fourth Bernoulli 
polynomial. For our purposes, a convenient explicit form is 



N - r 2 1 

7T 2 1 



k i ; / (^) = -t 4 E( 1 -^ 



AT 



_45 24tt 2 

where |A^- fe | 27r lies in the interval between and 2n. 
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\Ae jk \l(\A6 jk \ 27T -2nf 



(62) 



In many systems, broken symmetry phases are found at low temperatures and symmetry 
is restored at high temperatures. The phase structure of gauge theories as a function of 
temperature is unusual because the broken-symmetry phase is the high-temperature phase. 
A lattice construction of the effective action for Polyakov loops, valid for strong-coupling, 
is instructive [54TI5T] . The spatial link variables may be integrated out exactly if spatial 
plaquette interactions are neglected. Each spatial link variable then appears only in two 
adjacent temporal plaquettes, and may be integrated out exactly using the same techniques 
that are used in the Migdal-Kadanoff real-space renormalization group [HUES]- The resulting 
effective action has the form 



where K is a function of the lattice gauge coupling g 2 and the extent of the lattice in the 
Euclidean time direction n t , which is related to the temperature by n t a = 1/T. In the strong- 
coupling limit of the underlying gauge theory, the explicit form for K is K ~ (l/g 2 N) nt 
to leading order. In the weak-coupling limit, a Migdal-Kadanoff bond-moving argument 
gives K ~ 2N/g 2 n t . This effective action represents a Z(iV)-invariant nearest-neighbor 
interaction of a spin system where the Polyakov loops are the spins. It depends only on 
gauge-invariant quantities. Standard expansion techniques show that the Z(N) symmetry 
is unbroken for small K, and broken for K large. This model explains why the high- 
temperature phase of gauge theries is the symmetry-breaking phase: the relation between 
K and the underlying gauge theory parameters is such that K is small at low temperatures, 
and large at high temperatures, exactly the reverse of a classical spin system where the 
coupling is proportional to T~ x . For small values of nt, the deconffnement transition can be 
easily extracted, but the phase transition is in the strong-coupling region and far from the 
continuum limit. A systematic treatment of strong-coupling corrections has recently been 
shown to yield values for the critical lattice couplings (3 C = 2N/g 2 for 577(2) and SU(3) 
that are within a few percent of simulation results for 4 < N t < 16 [58J. For Abelian lattice 
gauge theories in (2 + l)-dimensions, a duality transformation maps the high-temperature 
symmetry-breaking phase of the gauge theory into the high-temperature unbroken phase of 
the dual spin system, and the low-temperature symmetry-breaking phase of the spin system 
into the unbroken phase of the gauge theory, consistent with the general behavior seen using 
strong-coupling arguments. 




(63) 



(Jk) 
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4- Universality 



The lattice construction of the Polyakov loop effective action is a concrete realization of 
Svetitsky-Yaffe universality |SH], which states that a second-order deconfinement transition 
in a (d + l)-dimensional gauge theory is in the universality class of classical spin systems 
in d dimensions with the same global symmetry. Lattice simulations indicate that all pure 
SU(N) gauge theories have a deconfining phase transition at some temperature Td, above 
which center symmetry is broken. In accordance with predictions based on universality, 
the deconfinement transition for an SU(2) gauge theory in 3 + 1 dimensions has been well- 
established as being in the universality class of the three-dimensional Ising model, exhibiting 
a second-order transition at Td- The deconfinement transition for SU(3) in 3 + 1 dimensions 
is first-order. This is consistent with Landau-Ginsburg predictions for a system with a Z(3) 
symmetry. The transitions for N > 3 appear to be first-order in 3+1 dimensions as well, with 
a smooth limit as N goes to infinity |60| |6T] . so the direct applicability of the universality 
argument in 3 + 1 dimensions is more limited than one might have expected. 

5. Equation of state 

Lattice simulations can be used to determine both the pressure p (T) and the internal 
energy e (T), along with the deconfinement transition temperature Td and related properties 
of the transition; see e.g. and references therein. For pure gauge theories below Td, both 
p and e are extremely small. This is expected because the mass of the lightest color singlet 
state, the scalar glueball, is substantially larger than Td- The pressure must be continuous 
at Td, but the first-order character of the deconfinement transition for N > 3 leads to a 
latent heat, indicated by a discontinuity in e. The thermodynamics quantities p and e as 
well as the order parameter (TrpP), show a rapid rise in the interval from Td to roughly 
3Td- Both p and e show a monotonic but slow rise towards their blackbody values at higher 
temperatures. 

The pressure p (T) and related thermodynamic quantities can be calculated in perturba- 
tion theory in terms of g 2 (T), the running coupling at temperature T. However, there is a 
barrier to perturbative calculations at 0(g 6 ), due to infrared divergences in the magnetic 
sector |63]. The best that can be done in perturbation theory is an expression for p(T) 
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that is valid up to O (g 6 log g) [64j . Comparison of different order of perturbation theory 
up to O (g 6 log g) shows a reasonable convergence only at temperatures greater than about 



6. Phenomenology of the deconfinement transition 

The availability of high-quality lattice data for the SU(3) pressure pjoTltjS] has led to a 
variety of attempts to model it phenomenologically. Deconfinement can be characterized 
broadly as a change in the number of degrees of freedom as the temperature is raised, 
and quark and gluon degrees of freedom manifest in thermodynamic behavior. The key to 
building a successful model is incorporating a mechanism for characterizing the change in 
degrees of freedom with temperature. Phenomenological models built around minimizing 
the free energy as a function of the order parameter are appealing and simple [4~2"j H?9"| 170] . 

One simple model |4"2| 171] follows from noting that the one-loop free energy for massive 
particles in a Polyakov loop background lead to expansions in which the first two terms have 
simple forms when M/T <C 1. Thus we arrive at a potential of the form 



It is important to note that the mass parameter M should not be interpreted as a gauge 
boson mass, nor do we limit ourselves to ML <C 1: the additional term in V g is purely 
phenomenological. The crucial feature of this potential is that for sufficiently large values 
of the dimensionless parameter M/T, the second term dominates and the potential leads to 
a Z(iV)-symmetric, confining minimum for P. On the other hand, for small values of M/T, 
the first term dominates and the pure gauge theory will be in the deconfined phase. This 
model has the nice property that Vg is a good representation of the gauge boson contribution 
for high temperatures. 

Minimizing Vq as a function of the eigenvalues of P leads to expressions for the pressure 
and associated thermodynamic quantities. In the case of SU(3), we can set the mass scale 
M by requiring that V g yields the correct deconfinement temperature for the pure gauge 
theory, with a value of Td = 270 MeV, giving M = 596 MeV. The pressure p is given by the 
value of —Vg at the minimum of the potential. The behavior obtained for the pressure, the 
energy density e, and the interaction measure A = (e — 3p) /T 4 are all roughly consistent 
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Figure 4: The Polyakov loop TrpP as a function of T in SU(3) using a phenomenological model 



with lattice simulations for T > Td and can be improved with the use of one or two additional 
free parameters |69 | I70 ] . As shown in Fig. |i| the order parameter TrP for the deconfmement 
transition jumps at Td, indicating a first-order deconfmement transition for SU(3). 

C. Role of quarks and other particles 

1. Chiral symmetry breaking 

If confinement is the most striking feature of QCD, chiral symmetry breaking is certainly 
a close second. If we imagine for a moment a variant of QCD in which all Nf quarks have 
the same mass, there will be a global SU (Nf) symmetry which rotates the different flavors 
into each other. In this theory, the left- and right-handed components of the Dirac fields, 
defined by 



i2l[7T]. 



qi/R = -(l ±7s)g 



(65) 



are connected by the Lagrangian mass term 



mqq = m (q L q R + q R q L ) 



(66) 
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which can be thought of as an interaction that flips the helicity of a quark. If the masses of 
the quarks are set to zero, then the left- and right-handed components are no longer coupled, 
and the symmetry group becomes SU(Nf)i £g> SU(Nf)n. At the classical level, there is also 
a U(1)a axial symmetry generated by 75, but it is anomalous. 

In our universe, the two lightest quarks, the u and d quarks, have masses m u and 
which are more than an order of magnitude smaller than Aqcd, so we are close to having an 
SU(2)l <g> SU(2)r symmetry. On the basis of models and, more recently, lattice simulations, 
we believe that in the massless limit this symmetry would be spontaneously broken to SU(2), 
the symmetry group of the massive theory, by the existence of a condensate (0\uu + dd\ 0). 
The breaking of this continuous global symmetry leads to the existence of three Goldstone 
bosons, the pions. 

Unfortunately, our ability to understand chiral symmetry breaking analytically has been 
limited. We assume that the basic mechanism is similar to the mechanism of supercon- 
ductivity. In the case of superconductivity, electron-hole interactions mediated by phonons 
lead to an expectation value for a composite field, while in gauge theories, the gauge bo- 
son interactions somehow give rise to some effective four-fermion interaction which leads 
to the formation of a qq condensate. However, there is no convincing derivation of this 
effective interaction from the gauge theory. In the absence of a more fundamental under- 
standing, Nambu-Jona Lasinio (NJL) models are often used to describe the chiral symmetry 
effects. In NJL models, a four-fermion interaction induces chiral symmetry breaking; the 
four-fermion interaction is non-renormalizable, and NJL models cannot be regarded as fun- 
damental. There has been a great deal of work on NJL models, both as phenomenological 
models for hadrons and as effective theories of QCD [72| 173] . 

NJL models use purely fermionic interactions as a proxy for the gauge theory interactions 
that give rise to chiral symmetry breaking. Typically, the relation between the gauge theory 
and an associated NJL model is fixed by matching important hadronic parameters such as 
f n . In the case of NJL models at finite temperature, and particularly PNJL models, it is 
common to assume that the NJL model parameters are fixed by T = hadronic parameters 
and remain constant as T is increased, at least up to the deconfinement temperature. 

As an example of an NJL model consider the representative model [71] 

L NJL = *p(ii-d-m )ij+^ [(VA» 2 + (^7 5 A>) 2 ] + g D [det^ (1 - 75) V + h.c] (67) 
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where ip is associated with Nf flavors of Dirac fermions in the fundamental or adjoint 
representation of the gauge group SU(N). The A a 's are the generators of the flavor symmetry 
group U(Nf); gs represents the strength of the four-fermion scalar-pseudoscalar coupling 
and go fixes the strength of an anomaly induced term. For simplicity, we can take the 
mass matrix mo to be diagonal: (mo)^ = mojSjk- If mo and gp are taken to be zero, 
Lnjl is invariant under the global symmetry U(Nf)i x U(Nf)n. With a nonzero, flavor- 
independent mass m and gjj ^ 0, the symmetry is reduced to SU(Nf)y x U(l)y. This 
model can be solved approximately in a self-consistent manner. The effective potential Vpo 
can be calculated at zero temperature as a function of the chiral condensates <7j = ijjji/jj : 

Nt 

V F0 K-) = £ 98*$ + 2g D {N f - 1) a,- - 2JV £ /" f^U? (68) 
j j i=i J { ' 

where the last term representing the sum of fermionic zero-point energies from the functional 
determinant. The zero-point energy of each mode uj^ = ^Jk 2 + m|, is written in terms of 
the constituent mass rrij, given in this model as 

m . = m o _ 2 g s a i - 2g D (N f - 1) JJ a h . (69) 

The sum over zero-point energies is divergent, and requires regularization. Differentiation 
with respect to a 3 - yields a set of gap equations which can be solved numerically for the 
chiral condesate. 



2. Quarks and deconfinement 

The addition of quarks or other particles to a gauge theory can completely change the 
theory's finite temperature behavior. Such particles can directly affect deconfinement, be- 
cause they can alter or even destroy center symmetry. These effects can be seen directly in 
the effective potential. 

The calculation of the one-loop effective potential for gauge bosons can be extended to 
other particles such as quarks in a very general way |74| . Because the effective potential 
is central for subsequent discussion, we treat the calculation in detail. All of the required 
expressions can be obtained from the one-loop effective potential for a charged scalar boson 
of mass M moving in a background U(l) Polyakov loop in d spatial dimensions. This is 
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given by 

V B (9) = \ [ — H In [1 - e -^ +ie ] + \ ( — d In [l - e'^ 9 ] (70) 

where tOk is given by \Jk 2 + M 2 . Now consider a scalar boson in the fundamental repre- 
sentation of an SU (N) gauge group moving in a uniform Polyakov loop background. The 
one-loop contribution to the finite-temperature effective potential is 

V = \ I (f^ ^ (1 " Pe '"" k) + ^ (1 " ^^"^ (71) 
where the first term is due to particles and the second term is due to antiparticles. A global 
unitary transformation puts P into the diagonal form 

Pjk = 6 jk exp(i9j) (72) 
and the partition function can be written as 

f = £M*;)- (73) 

3 

As a second example, consider the case of the gauge bosons themselves, which lie in the 
adjoint representation of the gauge group. The Polyakov loop in the adjoint representation 
is an (iV 2 — 1) x (iV 2 — 1) matrix. The partition function for the iV 2 — 1 particles is 

1 N 1 

j,k=l 

where the Sjk removes a singlet contribution, and the factor of 1/2 corrects for over-counting 
since Vb has both a particle and antiparticle contribution. The factor s accounts for spin 
degeneracy; in 3 + 1 dimensions s = 2, a consequence of the two possible polarization states 
of gauge bosons. 

For our third and final example, consider the evaluation of fermionic partition functions, 
which can be reduced to the general bosonic problem. A typical fermionic contribution of 
particle and antiparticle has the form 

V F (9) = -\ [ —. In [1 + e -^ +id ] - \ [ In [l + e^^l (75) 
PJ (27r) d L 1 P J (2iv) d 1 

which is easily written as 

V F (6) = -V B (tt + 9). (76) 
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For fermions in the fundamental representation of SU(N), the partition function is 

*£M*i) = + (77) 

where the factor s again accounts for spin degeneracy. 

A low-temperature expansion for Vb (9) can be generated for arbitrary spatial dimension 
d by expanding the logarithm and integrating term by term: 

M d/2+l/2 



V B {0) 



2d/2-3/2 7[ dj2+l/2pd/2+l/2 



E „d/2 + l/2 ^+l)/2 M M ) COS M 



(78) 



n=l 



which gives for scalar bosons in the fundamental representation 

M d/2+l/2 



V B {9) 



2d/2-l/2 7r d/2+l/2 / gd/2+l/2 n d/2+l/ 



m E ^m K <*w (^ M ) ( pn + ptn ) ( 79 ) 



with similar results for bosons in other representations. For fermions, we have 

M d/2+l/2 ^ (_-Qr 



2^/2-3/2^/2+1/2 ^d/2+1/2 



(80) 



n=l 



In a path integral representation, the factors of (—1)™ are a consequence of fermionic an- 
tiperiodic boundary conditions. 

A high-temperature expansion for Vb {9) may also be obtained in 3 + 1 dimensions 



V B {9) 
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(81) 
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M 4 
16^ 



In 



/3M 
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where # + is 9 modulo 2tt such that < 9 < 2tt. The sum ^ over I is over all integers with 



R G8M, e,0 = | [( + (^ - 2vr/) 2 ] 3/2 - 1 10 - 27T/I 3 - 1 1(9 - 2/vrl /3 2 M 2 - 
except that the divergent term is omitted when 1 = 0. 



(82) 



The first term in eqn. 81 is the blackbody free energy for two degrees of freedom, and 
depends only on the temperature and the angle 9. The second term, which is the leading 
correction due to the mass M at high temperatures, often appears in discussions of symmetry 
restoration at high temperatures with 9 = 0. The third term is closely associated with 
the n = Matsubara mode, which is the most infrared singular contribution to a finite 
temperature functional determinant. This term is responsible for non-analytic behavior in 



35 



finite temperature perturbation theory via the summation of ring diagrams. For example, 
in a scalar theory it gives rise to the A 3 / 2 contribution to the free energy; in QED, the 
contribution is e 3 [75] • The last term is logarithmic in the dimensionless combination /3M 
and independent of 9. In calculations of effective potentials, it typically combines with 
zero-temperature logarithms in such a way that the temperature T sets the scale of running 
coupling constants at high T . 

3. Lattice results 

Physical QCD has neither center symmetry, because there are quarks, not chiral sym- 
metry, because there is a quark mass term in the Lagrangian. Nevertheless, simulations of 
SU(3) gauge theories with quarks have revealed a complicated phase structure, where the 
number of flavors and the quark masses are crucial parameters. Quark masses span several 
orders of magnitude, from the light u and d quarks, with u, d -C ^qcd to the very heavy t 
quark with m t ^> Aqcd- Because m s is roughly on the order of A.qcd , lattice simulations 
aimed at realistic descriptions of hadronic physics are generally carried out with 2 + 1 flavors 
for sea quarks, that is, two light flavors and one intermediate relative to Aqcd- 

For gauge theories at finite temperature, the static approximation is the next logical step 
beyond pure gauge theories, and represents very well the effect of very heavy particles. In 
the case where quarks or other particles in the fundamental representation satisfy /3M ^> 1, 
we need keep only the leading term in the low-temperature expansion 



This generalizes immediately to particles in a representation R. Recalling the spin model 
interpretation, we see that a heavy particle behaves like an external field Kr coupled to the 
Polyakov loop TtrP in the representation R. For fermions in the fundamental representa- 
tion, this breaks Z(N) symmetry explicitly. In the case of SU(3) and other theories with 
first-order deconfinement transitions, the addition of a small Z(N) symmetry-breaking field 
does not remove the transition, as it would in the case of a second-order transition. Instead, 
a small symmetry-breaking field with hp > makes the deconfined phase more favorable 
at the deconfinement temperature, and there is a new deconfinement temperature, giving 
rise to a critical line in the T — hp plane. On both sides of the line, TrpP is non-zero, but 




(83) 
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there is a discontinuity that grows smaller with increasing hp. This behavior persists until 
a critical value of hp is reached, where the first-order line terminates in a critical end point, 
believed to be in the usual 4 universality class [751 - 178j . 

As quark masses are lowered, chiral symmetry effects become important. Phase tran- 
sitions, if present, lie at temperatures below the deconfinement transition temperature of 
the pure gauge theory, and the phase structure is very sensitive to quark masses. The re- 
sults of the extensive lattice studies that have been done is usually displayed in a so-called 
Columbia plot, which shows the phase structure of three-flavor QCD, with up, down and 
strange quarks, by showing lines of phase transitions as functions of a common up and down 
quark mass m u d and the strange quark mass m s . A version of this plot is shown in Figure 
|5j There are several interesting limiting cases contained in this diagram: 

• m ud — m s — oo: Because the quarks are completely removed from the dynamics, 
this limit is a pure gauge theory. For SU(3), the phase transition is first-order and 
associated with the spontaneous breaking of center symmetry, as described above. If 
the quark masses are decreased from infinity, center symmetry is explicitly broken, 
but the phase transition persists for large values of the quark masses, manifesting as a 
jump in (TrpP) from one non-zero value to another. As explained above, this region 
of first-order transitions terminates in what appears in the figure as a second order 
line. 

• r^ud — ; m s = oo This is the two-flavor chiral limit of QCD, with a second-order 
transition in the O (4) universality class that extends down from m s = oo until it 
meets another critical line at a tri critical point. 

• m ud — m s — This is the three-flavor chiral limit of QCD. Of course, it is less realistic 
than the two-flavor limit because the strange quark is sufficiently heavy to be on the 
order A. This is a first-order transition that extends into the m u d — m s plane before 
terminating in a second-order line. 

Physical QCD does not appear to have a finite temperature phase transition. Nevertheless, 
remnants of both the deconfinement transition of the pure gauge theory and the chiral 
transition of massless quarks remain. Near the critical temperature of the chiral transition 
in the massless theory, there is a crossover region. This is a relatively narrow range of 
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Figure 5: Schematic version of the Columbia plot showing phase structure as a function of the m u d 
and m s quark masses 

temperatures over which p, e, ipip and TtfP change rapidly. In physical QCD, this crossover 
marks the change from low-temperature hadronic behavior to high-temperature quark-gluon 
plasma behavior. 

4- connection of chiral symmetry breaking and confinement 

If a system has two or more order parameters, their behavior is typically correlated. This 
can be understood as arising from terms in the effective potential allowed by symmetry that 
couple the order parameters. In general, a first-order transition observed as a discontinuity 
in one order parameter will also give rise to a discontinuity in others. In the case of QCD 
and related theories, lattice simulations with fundamental representation fermions show that 
deconfinement and chiral symmetry restoration are strongly correlated. 
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The connection of chiral symmetry breaking and deconfmement can be seen in strong- 
coupling lattice models |79j . Integration over the spatial link variables does two things: 
generates the nearest-neighbor Polyakov loop coupling that drives the deconfmement tran- 
sition in strong coupling, and produces a four-fermion coupling that leads to spontaneous 
chiral symmetry breaking at low T. In order to explore the interrelationship of confine- 
ment and chiral symmetry breaking in a continuum model, we can use a generalization of 
Nambu-Jona Lasinio models known as Polyakov-Nambu-Jona Lasinio (PNJL) models [80J; 
see |81| for an alternative approach. NJL models have been used to study hadronic physics 
at finite temperature, but they include only chiral symmetry restoration, and do not model 
deconfmement. This omission is rectified by the PNJL models, which include both chiral 
restoration and deconfmement. In PNJL models, fermions with NJL couplings move in 
a nontrivial Polyakov loop background, and the effects of gluons at finite temperature is 
modeled in a semi-phenomenological way. 

The key idea of PNJL models is familiar: all particles move in a constant Polyakov loop 
background, and the effective potential will be a function of both the chiral order parameter 
ipip and P. The coupling between the chiral condensate and the Polyakov loop arises because 
the covariant derivative replaces the conventional derivative in the fermion kinetic 
term. However, contains only a constant gauge field background that gives rise to a 
non-trivial Polyakov loop. The complete PNJL effective potential consists of three terms 

Vp N j L = V G (P) + V F0 (aj) + V FT (P, Uj) (84) 

where Vq (P) reproduces the behavior of the pure gauge theory, Vp (°j) is the NJL effective 
potential at T = 0, and Vpp {P, &j) is the T^O part of the fermion effective potential, given 
by 

ln(l+Pe-^ )+h.c] (85) 



where the u)j~! s are defined in terms of the constituent masses as they were for the NJL 



model in Section III C If Vq is given by eqn. 64 and Vpo by 68 , then minimizing Vpnjl as 



a function of a for two flavors with tuq = 5.5 MeV gives the behavior shown in Fig. [6] \71\ . 
This shows the explanatory power of PNJL models. The constituent mass m is heavy at 
low temperatures, due to chiral symmetry breaking. The larger the constituent mass, the 
smaller the Z(3) breaking effect of the fermions, reflected in the small value of (TrpP) at 



39 



1.0 



m(T)/m(0) 




0.8 



0.6 



0.4 



0.2 



<TrP(T)>/3 



0.0 



50 



100 



150 
T (MeV) 



200 



250 



300 



Figure 6: Order parameters in a PNJL model of QCD with two light flavors as a function of 
temperature |71j . 

low temperatures. On the other hand, a small value for (TrpP) reduces the effectiveness of 
finite-temperature effects in restoring chiral symmetry. These synergistic effects combine in 
the case of fundamental representation fermions to give a single crossover temperature at 
which both order parameters are changing rapidly, consistent with the rapid crossover see 
in lattice simulations. 

Chiral and deconfmement behavior are not always as closely correlated as they are with 
quarks in the fundamental representation. In the case of adjoint fermions, lattice simulations 
have shown that chiral symmetry is restored at a substantially higher temperature than the 
deconfmement temperature, with the ratio of critical temperatures T c /Td around 7.8 [82, 83J. 
This behavior can also be reproduced in adjoint PNJL models |71| . 

D. finite density QCD 

The phase diagram of QCD in the temperature-chemical potential plane is of great in- 
terest. Unfortunately, lattice gauge theory is not directly applicable because of the sign 
problem. The nature of the problem can be seen from the one-loop effective potential. The 
effective potential Vft^ for a single flavor quarks at temperature T = and chemical 
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potential fi in a static Polyakov loop background P is 



/ 



d 3 k 



^Tr F [In (1 + Pe^-P" k ) + In (l + pt e -^-^fc)] 



(86) 




2T 



(27T) 



reflecting the fact that a non-zero chemical potential behaves as an imaginary U(l) Polyakov 
loop. This expression is complex for \i 7^ when P is non-trivial. This can be seen very 
clearly for heavy quarks, where Vpt> can be written as 



Because TrpP is complex for N > 3, the effective action for the gauge fields is complex. 
This is one version of the sign problem for gauge theories at finite density: the Euclidean 
path integral involve complex weights. This problem is a fundamental barrier to lattice 
simulations of QCD at finite density [84], 185] and a significant problem in other fields |86j . 

A large part of the interest in finite-density QCD lies in the possibility of a color supercon- 
ducting phase of hadronic matter in the interior of neutron stars j 16 j . At non-zero density, 
there are additional order parameters reflecting the possibility of Cooper pair formation 
whenever an attractive channel exists between fermions. As we have seen, chiral symmetry 
breaking can be viewed as the formation of a condensate of qq pairs in the vacuum. Other 
bifermion condensates are possible when the chemical potential is non-zero, giving rise to 
what is known as color superconductivity. This is an extension of the BCS mechanism from 
the formation of Cooper pairs of electrons to qq condensates. The mechanism for BCS su- 
perconductivity is an attractive force between Cooper pairs mediated by phonons. Early 
work on color superconductivity considered both interactions mediated by gluon exchange 
and induced by instantons via the axial U(l) anomaly. In practice, it is common to assume 
an NJL model with a set of four-fermion interactions, perhaps with an anomaly-induced 
term as well that is a 2Nj fermion coupling. 

Because quarks carry both color and flavor, it is necessary to consider both when dis- 
cussing appropriate order parameters for color superconductivity. Quarks carry flavor, color 
and spinor indices, so the quark condensate in a given channel can be characterized by 



where a and f3 are color indices, j and k are flavor indices, and a and b are spinor indices. The 
color-flavor-spin matrix characterizes a particular pairing, and the gap parameter A gives the 



Vft» ~f q ~ ~h F [e^Tr F P + e~^Tr F P + ] . 



(87) 




(88) 
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magnitude of the condensate in this channel. The condensate may be momentum-dependent. 
The familiar BCS condensate is a special case where the condensate is independent of posi- 
tion, so the pairing is between fermions of equal and opposite momentum, and a spin singlet. 
Unlike the fermion-antifermion condensates, the fermion-fermion condensates are not gauge 
invariant, and thus are similar to the Higgs scalar expectation value, which similarly depends 
on the choice of gauge. As with the Higgs mechanism, physical observables like the spectrum 
are gauge-invariant. There are strong arguments why the color-flavor locking phase (CFL) 
phase is favored for high-density QCD, combining antisymmetry in color and flavor with a 
Lorentz scalar behavior, giving 

(#C 75 <) oc e a ? A e ijA (89) 

where C is the charge-conjugation matrix acting on Dirac spinors. For a recent review of 
our current understanding of the phase diagram of dense QCD, see [87J. 



IV. PHASES OF GAUGE THEORIES ON R 3 x S* 1 



In the last few years, it has proven possible to construct four- dimensional gauge the- 
ories for which confinement may be reliably demonstrated using semiclassical methods 
[U |2j. These models combine Z(N) symmetry, the effective potential for P, instantons, 
and monopoles into a satisfying picture of confinement for a special class of models. All 
of the models in this class have one or more small compact directions. Models with an 
R 3 x S 1 topology have been most investigated, and discussion here will focus on this class. 
The use of one or more compact directions will cause the running coupling constant of an 
asymptotically free gauge theory to be small, so that semiclassical methods are reliable. For 
example, if the circumference L of S 1 is small, i.e., L <C A -1 , then g(L) <C 1. However, this 
leads to an immediate problem: generally speaking, one or more small compact directions 



lead to breaking of Z(N) symmetry in those directions. As we have seen in Section III B 
for the case of finite temperature gauge theories, where L = (3 = 1/T, the effective potential 
for the Polyakov loop is easily calculated to lowest order in perturbation theory; for a pure 



gauge theory it is given by eq. ( 56 ) 



7T^ 

n=l 
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where the trace of P in the adjoint representation is given by Tr^P = \TrpP\ 2 — 1. This 
effective potential is minimized when P = zl where z G Z(N), indicating that Z(N) 
symmetry is spontaneously broken at high temperatures where the one-loop expression is 
valid. However, it is possible to maintain Z(N) symmetry even when L is small by modifying 
the action. This leads to a perturbative calculation of possible phase structures, which turns 
out to be very rich, as well as a perturbative understanding of Polyakov loop physics in 
the confined phase. Furthermore, the restoration of Z(N) symmetry for small L leads to 
a non-perturbative mechanism for confinement, as measured by Wilson loops orthogonal 
to the compact direction. In this confinement mechanism, a key role is played by finite- 
temperature instantons, also known as calorons, and their monopole constituents. Thus we 
obtain a realization of a long-held scenario for quark confinement, based on ideas originally 
proposed by Mandelstam jHH [89] and 't Hooft [ffl EE]. 



A. Restoring center symmetry via deformations or periodic adjoint fermions 

There are two broad approaches to maintaining Z(N) symmetry for small L. The first 
approach deforms the pure gauge theory by adding additional terms to the gauge action 
[21 [H2J US]- The general form for such a deformation is 



/oo 
d 3 x a kTr A P i 



x,Xi) k (91) 



where the value of x^ is arbitrary and can be taken to be 0. Such terms are often referred 
to as double-trace deformations; see Fig. [7j If the coefficients at are sufficiently large, 
they will counteract the effects of the one-loop effective potential, and Z(N) symmetry 
will hold for small L. Strictly speaking, only the first [N/2] terms are necessary to ensure 



confinement. As discussed in Section [HI B[ it is easy to prove that for a classical Polyakov 
loop P, the conditions TrpP k = with 1 < k < [N/2] determine the unique set of Polyakov 
loop eigenvalues that constitute a confining solution, i.e., one for which TtrP = for all 
representations with kn ^ |42j . The explicit solution is simple: up to a factor necessary to 
ensure det P = 1, the eigenvalues of P are given by the set of iV'th roots of unity, which are 
permuted by a global Z(N) symmetry transformation. The effective potential associated 
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Figure 7: A double-trace Polyakov loop on ii 3 x S 1 , 
with S is given approximately by 

— 9 TV , , P n _ 2 : 

K// (P, /5) = zJiE + £ (92) 

^ n=l fe=l 

For sufficiently large and positive values of the a^'s, the confined phase yields the lowest 
value of V e ff. However, a rich phase structure emerges from the minimization of V e ff for 
intermediate values of the coefficients a*.. For N > 3, the effective potential predicts that 
one or more phases may separate the deconfined phase from the confined phase. In the 
case of SU(3), a single new phase is predicted, and has been observed in lattice simulations 
|2j. For larger values of N, there is a rich set of possible phases, including some where 
Z(N) breaks down to a proper subgroup Z(p). In such phases, particles in the fundamental 
representation are confined, but bound states of p particles are not [92J. 

Lattice simulations of SU(3) and SU{4) agree for small L with the theoretical predictions 
based on effective potential arguments [2j. The phase diagram of SU(3) as a function of 
T = L~ l and a\ has three phases: the confined phase, the deconfined phase, and a new 
phase, the skewed phase. In general, the three phases of the eigenvalues of the Polykov 
loop may be taken to be the set {^1,^2,^3} where 6\ + 6 2 + #3=0. For all three phases, 
it is possible to use Z(3) symmetry to make TrpP real, and reduce the phases to the set 
{0,6,-6} such that TrpP = 1 + 2cos#. The deconfined phase is represented by 6 = 0, 
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the confined phase is given by 9 = 2tt/3, and the skewed phase by 6 = tt. An important 
result obtained from the lattice simulation of SU(3) is that the small- L confining region, 
where semiclassical methods yield confinement, are smoothly connected to the conventional 
large- L confining region. In the case of 577(4), a sufficiently large value of a\ leads to a 
partially-confining phase where Z(A) is spontaneously broken to Z(2). Particles with k = 1 
are confined in this phase, i.e., (TtfP {x)) = 0, but particles with k = 2 are not, as indicated 
by (TrpP 2 (x)) ^ 0. As in the case of 577(3), TrpP can be made real. Perturbation theory 
then predicts a deconfined phase where the phases of the eigenvalues of P are {0, , 0, 0, 0} , a 
confined phase where they are {vr/4, 37r/4, 57r/4, 77r/4}, and a partially confined phase where 
the phases are {7r/2,7r/2, — n/2, —it/2}. The confined phase corresponds to the matrix Pi 



and the partially confined phase to the matrix P2 described in Section III B 

Another approach to preserving Z(N) symmetry for small L uses fermions in the adjoint 
representation with periodic boundary conditions in the compact direction [T|. In this case, 
it would be somewhat misleading to use /3 as a synonym for L, because the transfer matrix for 
evolution in the compact direction is not positive-definite. Periodic boundary conditions in 
the compact direction imply that the generating function of the ensemble, i.e., the partition 
function, is given by 



Tr 



■l) F e- LH 



(93) 



where F is the fermion number and H is the Hamiltonian in the compact direction. 
This graded ensemble, familiar from supersymmetry, can be obtained from an ensemble 
Tr [exp (/3/j,F — 0H)] with chemical potential \x by the replacement /3fi — > in. This system 
can be viewed as a gauge theory with periodic boundary conditions in one compact spatial 
direction of length L = /3, and the transfer matrix in the time direction is positive-definite, 
The use of periodic boundary conditions for the adjoint fermions dramatically changes 
their contribution to the Polyakov loop effective potential. In perturbation theory, the 
replacement — > in shifts the Matsubara frequencies from 0u n = (2n + 1) 7T to (3u n = 2nir. 
The one loop effective potential is now essentially that of a bosonic field, but with an overall 
negative sign due to fermi statistics [71]. The sum of the effective potential for the fermions 
plus that of the gauge bosons gives 

K // (P,/?,™,A' / ) = ^E l!: f f! 

H n=l 

where Nf is the number of adjoint Dirac fermions and m is their mass [49J. Note that the 



2N f (3 2 m 2 K 2 (n/3m) 



ri 2 



(94) 
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first term in brackets, due to the fermions, is positive for every value of n, while the second 
term, due to the gauge bosons, is negative. 

The largest contribution to the effective potential at high temperatures is typically from 
the n = 1 term, which can be written simply as 

[2N f (3 2 m 2 K 2 (/3m) - 2] [\Tr F P\ 2 - l] (95) 

where the overall sign depends only on Nf and (3m. If Nf > 1 and (3m is sufficiently small, 
this term will favor TrpP = 0. On the other hand, if /3m is sufficiently large, a value of P 
from the center, Z(N), is preferred. Note that an M = 1 super Yang-Mills theory would 
correspond to Nf = 1/2 and m — 0, giving a vanishing perturbative contribution for all n 
[9"4"t 193] . In that case, non-perturbative effects lead to a confining effective potential for all 
values of /3. In the case of Nf > 1, each term in the effective potential will change sign in 
succession as m is lowered towards zero. For larger values of N, this leads to a cascade of 
phases separating the confined and deconfined phases |48| . Numerical investigation shows 
that the confined phase is obtained if N /3m < 4.00398 [49]. As m increases, it becomes 
favorable that TrpP n ^ for successive values of n. If N is even, the first phase after the 
confined phase will be a phase with Z(N/2) symmetry. As m increases, the last phase before 
reaching the deconfined phase will have Z(2) symmetry, in which k = 1 states are confined, 
but all states with higher k are not. Lattice simulations of SU(3) with periodic adjoint 
fermions are completely consistent with the picture |96j predicted by the effective potential, 
with a skewed phase separating the confined phase and deconfined phase. For N > 3, there 
are generally phases intermediate between the confined and deconfined phases which are 
not of the partially-confined type. Careful numerical analysis appears to be necessary on a 
case- by-case basis to determine the phase structure for each value of N |48j . There has been 
very interesting work on lattice gauge theories in the large- N limit that indicates the same 



rich phase structure [971 198] ; this work will be discussed in Section VI 



There are some interesting additional issues arising when periodic adjoint fermions are 
used to obtain Z(N) symmetry for small L. There are strong indications from strong- 
coupling lattice calculations jTS] and the closely related Polyakov-Nambu-Jona Lasinio 
(PNJL) models |SU] that the mass m that appears in the effective potential V(P) should 
be regarded as a constituent mass that includes the substantial effects of chiral symmetry 
breaking. In strong-coupling lattice calculations and PNJL models, this effect is responsible 
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for the coupling of P and i/ji/j in simulations of QCD at finite temperature. Thus it is impor- 
tant that simulations of SU (3) with Nf = 2 flavors of adjoint fermions show explicitly that 
the Z(iV)-invariant confined phase is regained when the fermion mass is sufficiently small 
|96| . However, this raises another issue. For a simple double-trace deformation, the small- L 
and large-!/ regions are smoothly connected. A semi-phenomenological analysis based on 
a PNJL model [71j suggests that this is also the case with periodic adjoint fermions, but 
additional modifications of the action may be necessary to realize the connection. 

B. Monopoles and instantons on R 3 x S 1 

The non-perturbative dynamics of confining gauge theories on R 3 x S 1 are based on 
Polyakov's analysis of the Georgi-Glashow model in three dimensions [HH]- This is an SU(2) 
gauge model coupled to an adjoint Higgs scalar. The models we are considering thus differ by 
the addition of a fourth compact dimension and a change to the action designed to maintain 
Z(N) symmetry. The four-dimensional Georgi-Glashow model is the standard example of 
a gauge theory with classical monopole solutions when the Higgs expectation value is non- 
zero. These monopoles make a non-perturbative contribution to the partition function Z. 
In three dimensions, these monopoles are instantons. Polyakov showed that a gas of such 
three-dimensional monopoles gives rise to non-perturbative confinement in three dimensions, 
even though the theory appears to be in a Higgs phase perturbatively. 

Because L is small in the R 3 x S l models we consider, the three-dimensional effective 
theory describing the behavior of Wilson loops in the non-compact directions will have 
many features in common with the three-dimensional theory first discussed by Polyakov. 
In the four-dimensional theory, monopole solutions with short worldline trajectories in the 
compact direction exist, and behave as three-dimensional instantons in the effective theory; 
see Fig. [8] In models on R 3 x S 1 , the role of the three-dimensional scalar field is played 
by the fourth component of the gauge field A4. In a gauge where the Polyakov loop is 
diagonal and independent of £4, P has a vacuum expected value induced by the perturbative 
effective potential. However, there is another way to understand the presence of monopoles 
in this phase, based on studies of instantons in pure gauge theories at finite temperature 
and the properties of the KvBLL caloron solution [100- 102]. If the Polyakov loop has a 
non-trivial expectation value, finite-temperature instantons in SU(N) may be decomposed 
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Figure 8: Short monopole worldline on i? 3 x S . 



into N monopoles, and the locations of the monopoles become parameters of the moduli 
space of the instanton. In the case of SU(2), an instanton may be decomposed into a 
conventional BPS monopole and a so-called KK (Kaluza-Klein) monopole. The presence of 
the KK monopole solution differentiates the case of a gauge field at finite temperature from 
the case of an adjoint scalar breaking SU(N) to [/(l)^" 1 , in which case there are N — 1 
fundamental monopoles. We will consider in detail the simplest case of N = 2. 

The BPS monopole is found using the using the standard arguments [103] 1104] . The 
Euclidean Lagrangian £ can be taken to be 

C = \ {F^f + +V eff (P) (96) 

where V e ff includes both the one- loop gluonic effective potential and the additional term 
that prevents Z(N) symmetry breaking. This can also be written as 

C = \ (DjAi) 2 + \ (B,) 2 + V eff (P) . (97) 



We can associate with C an energy defined by 



E= d 3 x 



\m 2 + \(D 3 A A f + Veff{P) 



(98) 



as well as an action S = LE. We will concern ourselves for now with the solutions in the 
BPS limit, in which the effective potential V e ff is neglected, but the boundary condition on 
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P at infinity imposed by the potential is retained. We can write the energy as 



E = (fx 



(99) 



This expression is a sum of squares plus a term which can be converted to a surface integral, 
giving rise to the BPS inequality 

E > =f J dSjBjAi. (100) 

The BPS inequality is saturated if the equality Bj = ^fDjA^ holds. For the case of a single 
monopole at the origin, we require the fields at spatial infinity to behave as 

lim A\ = w — 

r— >oo r 

lim A" = e ai] ^. (101) 



r—toc ' gr 2 



Note that w is related to the eigenvalues of P at large distances by w = 29/ gL. Note that 
A4 has the usual hedgehog form. A? is chosen such that covariant terms vanish at infinity: 
(DiA^" 1 = 0. With the 't Hooft-Polyakov ansatz, the general expressions for the fields 
become 

A\ = wh (r) — 
r 

gr 2 

where we define w > and require h(oo) = 1 or — 1, and a(oo) = 1 to obtain the correct 

asymptotic behavior. We must also have /i = a = 0atr = 0to have well-defined functions 

at the origin. We identify a magnetic flux 

f r a Ait 

$ = ± / dSjB?— = T— (103) 
J 3 r g 

where the + sign corresponds to the case h(oo) = 1 and — corresponds to h(oo) = — 1. The 
energy of the BPS monopole can be written as 

4.71W 

Ebps = T®w = . (104) 

9 

In addition to the BPS monopole, there is another, topologically distinct monopole which 
occurs at finite temperature when A4 is treated as a Higgs field [94j. Starting from a static 
monopole solution where \A±\ — w at spatial infinity, we apply a special gauge transformation 



U special exp 



171X4 q 
T 



(105) 
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where r l is the Pauli matrix. U spec i a i transforms in such a way that the value of A4 
at spatial infinity is shifted: w — > w — 2n/gL. If we instead start from a static monopole 
solution such that A4 = 2n/gL — w at spatial infinity, then the action of U specia i gives a 
monopole solution with A^ = —w at spatial infinity. A final constant gauge transformation 
U const = exp [i7TT 2 /2] yields a new monopole solution with A4 = w at spatial infinity. The 
distinction between the BPS solution, which is independent of x 4 , and the KK solution is 
made clear by consideration of the topological charge. The action of U specia i followed by 
U const increases the topological charge by 1 and changes the sign of the monopole charge. 
Thus the KK solution is topologically distinct from the BPS solution because it carries 



instanton number 1. The BPS antimonopole has magnetic charge opposite to the BPS 



monopole, and hence the same as that of the KK monopole. The KK monopole has the 
same magnetic charge as the BPS monopole, but carries instanton number —1. This is 
all completely consistent with the KvBLL decomposition of instantons in the pure gauge 
theory with non-trivial Polyakov loop behavior, where SU (2) instantons can be decomposed 
into a BPS monopole and a KK monopole. Our picture of the confined phase is one where 
instantons and anti-instantons have "melted" into their constituent monopoles and anti- 
monopoles, which effectively forms a three-dimensional gas of magnetic monopoles. In the 
BPS limit, both the magnetic and scalar interactions are long-ranged; this behavior appears 
prominently, for example, in the construction of 7V-monopole solutions in the BPS limit. 
The BPS solution has action 

S BPS = ^ = S 4. (106) 

9 g 2 

For the KK solution, we have instead 

4tt (2tt - gLw) 4tt (2tt - 29) 

OKK = K = n • l 1U 'J 

9 9 

The sum Sbps + Skk is exactly 8n 2 /g 2 , the action of an instanton. For 9 = n/2, the Z(2)- 
symmetric value for SU(2), Sbps — Skk- This extends to SU(N), where the action of a 
monopole of any type is 8n 2 /g 2 N. 

Although we used the BPS construction to exhibit the existence and some properties of 
the monopole solutions of our system, we must move away from the BPS limit to ensure 
that magnetic interaction dominate at large distances, i.e., that the three-dimensional scalar 
interactions associated with A4 and are not long-ranged. This behavior is natural in the 
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confined, where the characteristic scale of the Debye (electric) screening mass associated 
with A 4 is large, on the order of g/L. It is well known that the BPS bound for the monopole 
mass holds as an equality only when the scalar potential is taken to zero. Numerical studies 
[105J have shown that the monopole action is given in general for SU (2) as 

LE BPS C(e) (108) 

where C a function of the quartic term in the potential that varies from C = 1 in the 
BPS limit to a maximum value C (oo) = 1.787. Thus corrections to the BPS result for the 
monopole mass and action due to the potential terms are less than a factor of two. We will 
henceforth use the exact results for the actions in the BPS limit, neglecting corrections from 
V e ff for the sake of simplicity of notation. 

The 577(2) construction of BPS and KK monopoles extends to SU (N) in the standard 
way, via the embedding of SU(2) subgroups in SU(N). There are N — 1 BPS monopoles 
and 1 KK monopole inside an instanton. In the confined phase, each of the N monopoles 
has action 87f 2 /g 2 N. It has long been thought that instanton effects must be suppressed in 
the large- N limit, because instanton effects would vanish as exp (—cN) in the limit N — > oo 
with A = g 2 N fixed |106j . In contrast, we see that the effects of monopole constituents of 
instantons are not suppressed by the large- N limit. 

C. Monopoles and space-like string tensions 

It is important to understand that Wilson loops in planes orthogonal to the compact 
direction should show area law behavior, even if the Z(N) symmetry associated with the 
compact direction is broken. This is an old observation about the deconfined phase |107l 108J 
which is very clearly observed in lattice simulations of SU(2) and SU (3) at temperatures 
above the deconfinement transition |109l fllOj. At first sight, this seems to directly conflict 
with the association of deconfinement with the loss of area-law behavior for Wilson loops. 
However, the introduction of a compact direction, as in the case of finite temperature, explic- 
itly breaks space-time symmetry. In the case of finite temperature, Wilson loops measuring 
electric flux have perimeter behavior in the deconfined phase; Wilson loops measuring mag- 
netic flux still obey an area law. This asymmetry in behavior can be understood on the basis 
of center symmetry. The full center symmetry of an SU(N) gauge theory on a <i-dimensional 
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Figure 9: Spatial Wilson loop in R 3 x S 1 geometry. 

hypertorus T d is Z(N) d . While the Z(N) symmetry may break spontaneously in the short 
compact direction, the other Z(N) symmetries are unbroken, and thus the associated Wilson 
loops obey an area law. 

In order to understand the effects of monopoles play in the confined phase, we must 
analyze their interactions. We begin with a discussion of quantum fluctuations around the 
monopole solutions. The contribution to the partition function of a single BPS monopole 
at finite temperature was considered by Zarembo jlllj - The measure factor dfi a associated 
with the collective coordinates (moduli) of the monopole solution, including the Jacobians 
from the zero modes is given by |95| 



djjL a = /i 4 



d 3 



x 



- J r 



2tt 



;Jrh 



(109) 



(2nf 2 X Jo (2tt) 1/2 ' 

where x is the position and <ft the £7(1) phase of the monopole and fi is a Pauli-Villars 
regulator. The label a denotes the type of monopole, a = {BPS, KK, BPS, KK} . The 
Jacobians are 

J x = S 3 J 2 , = NLSl' 2 . (110) 

Each of the four zero modes contributes a factor of fj,. In the BPS limit, each monopole 
carries an overall factor 

Z a = cfi 7 / 2 (NL) 1/2 S 2 a exp[-S a + 0(l)] J d 3 x 

= £aexp[-S a ] j d 3 x (111) 
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in its contribution to Z . The factor £ a is cfi 7 ' 2 {NL) l/2 S 2 a where c is a numerical 
constant and the factor of d 3 x represents the integration over the location of the monopole. 
From the construction of the KK monopole, we see that we have ^kk (#) = £,bps (tt — 6). 

The renormalization of the functional determinant arising from quantum fluctuations 
around the monopole solution is particularly simple in the confined phase, as first observed 
by Davies et al. in the corresponding supersymmetric model [94]. The dependence on 
the Pauli-Villars regulator is removed, as usual, by coupling constant renormalization. The 
relation at one loop of the bare coupling and the regulator mass \x to a renormalization-group 
invariant scale A is 

where fe is the first coefficient of the /3 function divided by iV: 

11 4 n f C(R f ) 1 n b C(R b ) 

bo ~ y " 3 ' ~~iv e ' ^v - (113) 

where n/ is the number of flavors of Dirac fermions in a representation Rf, n b is the number 
of flavors of real scalars in a representation Rf,, and C(R) is obtained from Ttr (T a T 6 ) = 
C(R)5 ab . For the case of a pure gauge theory with a deformation, there are four collective 
coordinates and this gives a factor of /i 4 . The functional integral over gauge degrees of 
freedom gives rise to a factor det' [— D 2 ] 1 oc /i^ 1 / 3 and the action contributes a factor 
exp (— 87T 2 /g 2 N) in the confined phase. Thus the contribution of a single monopole to the 
partition function gives a factor 

^4-i e -8*WiV = ^H/Sg-S^iV = A ll/3_ (n4) 

Thus detailed calculation confirms what we might have guessed on dimensional grounds: 
the contribution £ a e~ 871-2 L~ 3 (AL) 11 / 3 . Note that the eliminations of renormalization- 

dependent quantities by renormalization-independent quantities depends crucially on the 
coefficient of 1/g 2 in the action. 

The interaction of the monopoles is essentially the one described by Polyakov in his 
original treatment of the Georgi-Glashow model in three dimensions [99] , slightly generalized 
to include both the BPS and KK monopoles. Let us consider, say, a BPS-type monopole 
and KK-type monopole located at x\ and X2 in the non-compact directions, with static 
worldlines in the compact direction. The interaction energy due to magnetic charge of such 
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a pair is 



4tt\ 2 1 



Ebps-kk = -[ — \ -tj^ ITT (115) 

V g J 47r \%i - x 2 \ 

and the associated action is approximately Sbps + Skk + LEbps-kk- As discussed above, 
this will be larger than the value obtained from the Bogomolny bound, but of the same order 
of magnitude. There is an elegant way to capture the dynamics of the monopole plasma, 
using an Abelian scalar field a dual to the magnetic field. Assuming that the Abelian 
magnetic gauge field is three-dimensional for small L, we may write 

d"x l -Bl = J d'x^^id.a) 2 (116) 

where the normalization of a is chosen to simplify the form of the interaction terms. The 
three-dimensional effective action is given by 



Leff = {d s <rf -J2t a e- S « +i ^ (117) 



f_ 
32vr 2 L 

a 

where the sum is over the set {BPS,KK,BPS,KK}. Each species of monopole has its 
own magnetic charge sign q a = ± as well as its own action S a . The coefficients £ a represent 
the functional determinant associated with each kind of monopole, but the combination 
£ a exp(— S a ) may be usefully regarded as a monopole activity in terms of the statistical 
mechanics of a gas of magnetic charges. The generating functional 



Z a = J [da] exp — J d 3 xL e ff 



(118) 



is precisely equivalent to the generating function of the monopole gas. This equivalence may 
be proved by expanding Z a in a power series in the £ a 's, and doing the functional integral 
over a for each term of the expansion. 

The magnetic monopole plasma leads to confinement in three dimensions. For our effec- 
tive three-dimensional theory, any Wilson loop in a hyperplane of fixed x±, for example a 
Wilson loop in the x\ — x 2 plane, will show an area law; see figure |9j The original procedure 
of Polyakov [99J may be used to calculate the string tension, where the presence of a large 
planar Wilson loop causes the dual field o to have a discontinuity on the surface associated 
with the loop and a half-kink profile on both sides. However, an alternative procedure is 
simpler where the discontinuity in the gauge field strength induced by the Wilson loop is 
moved to infinity so that the string tension is obtained from the kink solution connecting 
the two vacua of the dual field a [93J . 
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In the confined phase, the action and functional determinant factors for all four types 
of monopoles are the same, so we denote them by Sm and £m- The potential term in the 
mixed and confined phases then reduces to 

- J2 iae~ Sa+iqaU -> 4£ M e~ SM [1 - cos (a)} (119) 

a 

which has minima at a = and o = 2tt; we have added a constant for convenience such 
that the potential is positive everywhere and zero at the minima. A one- dimensional soliton 
solution a s (z) connects the two vacua, and the string tension a^d for Wilson loops in the 
three non-compact directions is given by 

/+00 
dzL eff (a z (z)) (120) 
-00 

which can be calculated via a Bogomolny inequality to be 



^= 4 {ffe^. (121) 

Of course, this result depends on L and cannot be used outside the region AL <C 1; never- 
theless, this is a concrete realization of confinement in a four- dimensional field theory via 
non-Abelian monopoles. 



D. The special case of supersymmetry 

The earliest indication that confinement on R 3 x S 1 was interesting and tractable came 
about in an unlikely way, in an effort to resolve a controversy over adjoint fermion con- 
densates in four- dimensional M = 1 supersymmetric SU(N) gauge theories. Davies et al. 
resolved the controversy by considering supersymmetric SU(N) gauge theories on R 3 x S 1 
|94 [ l95 ] . In such theories, supersymmetry requires that the number of adjoint fermion degrees 
of freedom ("gluinos") equals the number of gauge degrees of freedom, and that the adjoint 
fermions have the same boundary conditions as the gauge fields. This model interpolates 
between four-dimensional Af = 1 supersymmetric SU(N) gauge theories as L — > 00 and 
three-dimensional M = 2 supersymmetric SU(N) gauge theories in the limit L — > with no 
intervening phase transition as L is varied. In the three-dimensional limit, the scalar field of 
the d = 3 M = 2 theory is simply the n = Matsubara mode of A 3 . We follow closely the 
original notation and treatment of SU(2) by Davies et al. j9l], but see also [H I112j ; some 
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subtleties of the calculation that were previously unappreciated are discussed in detail in 
[113J. As before, a non- vanishing expectation value for A\ breaks the SU(2) symmetry to 
U (1), and the low-energy effective theory will have only £7(1) symmetry. The superpotential 
of this model is trivial in perturbation theory and the effective potential is completely non- 
perturbative. The superpotential can be written in terms of chiral and anti-chiral M = 1 
superfields 

$ = Z + V20 a V a + 9 a 9 a J : 

$ = z + V2d a v a + e a e a F 



where Z is a complex scalar <ft + i^j that includes both the three-dimensional scalar <fi = A\ 
and a field 7 which is dual to the £7(1) magnetic field £?| and thus proportional to o in our 
notation. The effective action is 
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where W (<&) is given by 
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The parameter M is M = 167r 2 L 2 Mp y with Mpy the Pauli-Villars regulator associated with 
the functional determinant. The effective potential is given by 
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(125) 

r / \ F F J\ 

The minimum of V e ff preserves supersymmetry: 7 ((Z) = 0) implies (0) = tt//3 and (7) = 0, 

and leads to TrpP = . The effective potential for 7 is then given by 



V,„ (</> = */ 
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(126) 



The string tension is again obtained from the one-dimensional kink solution of the equation 
of motion obtained from the effective action for the dual field. This model represents a 
lower boundary of where periodic adjoint fermions work to restore Z(N) symmetry, with 



56 



the superpartner of the gauge boson corresponding to Nf = 1/2 for the number of Dirac 
fermions. For Nf = 1/2, the one-loop effective potential Eqn. 94 favors the deconfmed 
phase for when the fermion mass is positive, and vanishes identically when it is zero. 



E. Polyakov loop string tensions 

Polyakov loop string tensions are calculable perturbatively in the small- L confining region 
from small fluctuations about the confining minimum of the effective potential |4"9l 1114] ; see 
also [93j. The effective Lagrangian can be written in terms of the phases of the Polyakov 
loop eigenvalues 9j 

1 - 

^E^f + ^^l (127) 
9 i=i 

where the kinetic term is obtained from the standard kinetic term for gauge theories and the 
potential term is the one-loop effective potential. In the high-T (small-L) confining region, 
the minimum of the effective potential will have Z(N) symmetry. A convenient form for 
such a solution is 

(P ) jk = wz^5 jk (128) 

where z = exp (27ri/N) and w is a phase factor that ensures detPo — 1- For small fluctua- 
tions, we write 

P = P e m (129) 

where 59 lies in the Cartan algebra of SU(N), i.e., is a diagonal traceless matrix, and 
therefore has N — 1 independent components. The moments of the Polyakov loop are 

N 

Tr F P k = w k z jk e ikse * (130) 
If k is not divisible by N, we have approximately 

N 

Tr F P k ~ ikw k z ik Mj = ikw k 56 k (131) 

3=1 

where 59 is the discrete Fourier transform of 58, related by 

N 

56 k = ^z jk 56 3 (132) 
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so^^Y^z-nh. (133) 

k=l 

Note that the reality of 58 implies that 88* k = 50N-k', the last Fourier component, 54>n, is 
identically zero, due the tracelessness of 5<f). If k is divisible by N we have instead 

N N 

Tr F P k ~N--k 2 J2 (S9j) 2 = N — — k 2 ]T {59 m 88 N -rn) ■ (134) 

j=l m=l 

These formulae allow us to write the three-dimensional effective action in terms of the 58k 
to quadratic order. For each value of k, the terms in the potential which contribute have 
n = k mod N, n = N — k mod N, or n = mod N. We obtain a different mass cr^/T 
for each Fourier component 56^ measured by using Polyakov loops of iV-ality k. There 
are [N/2] different string tensions because a k = <J^_ k . The string tensions obtained, for 
example, in lattice simulations will depend on the operators used. For example, the operators 
{TrpP) k ~ [f^y an< ^ Trp(P k ) ~ 58~k both have iV-ality k, but the two-point correlation 
function of the first operator will decay as exp [—kafr/T] , while the second operator decays 
as exp [— cr^r/T]. Group characters will in general have a complicated pattern of mixing. 
Consider the symmetric/anti-symmetric representations of dimension iV(iV ± l)/2, with 
characters given by 

Xs /a(P) = (1/2) [{Tr F Pf ± Tr F P 2 ] . (135) 

These operators will have two different contributions, corresponding to 2crf' ) and ■ At 
higher orders in perturbation theory, there will be mixing, and ultimately only the lightest 
state for a given iV-ality will be seen at large distances. This phenomena, often referred to 
as string breaking, was first observed in lattice simulations in SU(2) Higgs models |115H117] 
and later seen in pure SU(2) gauge theory in (2 + 1) dimensions [118J. 

The behavior seem in the models considered in this section is very different from the 
behavior seen in exactly solvable two-dimensional gauge theories or in lattice simulations 
of four-dimensional gauge theories. This is not surprising. For the case of double-trace 
deformations, there are at least [N/2] independent parameters which can be adjusted, and 
which can vary in turn the aj?. For the case of periodic adjoint fermion, the results are 
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continuous functions of /3m. The string tensions are of order g: 



T 



9 2 N 



Nftn 



V [K 2 ((k + jN)/3m) + K 2 ((N-k + jN)Pm) - 2K 2 ((j + l)Npm)\ 

IT* Z — / 

j=0 



-g 2 N 



-2 

T 2 
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irk 
TV 



- 1 



(136) 



where the gluon contribution has been summed in the last term. Note that the symmetry 



a k = a N-k i s manifest in this formula. The m = limit has the simple form 



(2N f -l)g 2 T 2 



3 csc^ 



irk 
~N 



1 



(137) 



T / 3N 

and is a good approximation for j3m 1. This scaling law is not at all like either Casimir or 
sine-law scaling, because the usual hierarchy a£ + 1 > a£ is here reversed. Because we expect 
on the basis of SU(3) simulations that the high-temperature confining region is continuously 
connected to the conventional low-temperature region, there must be an inversion of the 
string tension hierarchy between the two regions for all N > 4. Other features of the string 
tension behavior are in line with our expectations. For the case Nj = 1/2 , corresponding 
to a single multiplet of adjoint Majorana fermions, the perturbative string tension vanishes, 
and it is the non-perturbative contribution to the effective potential induced by monopoles 
that gives rise to the string tension in this case |9"H [95] . The large- iV limit of Eqn. 137 is 
smooth. For fixed k as iV — > oo, we have 



T 



(2N f - 1) AT 2 

7T 2 k 2 



(138) 



where A is the 't Hooft coupling g 2 N . The Polyakov loop string tensions obtained in the 
small-!/ confining region, while conforming to general principles, do not appear to tell us 
much about confinement at large L. 



F. Other geometries 

It is possible to apply the same methods used for gauge theories on R 3 x S* 1 to other 
geometries, such as S 1 x S 3 or R 2 x T 2 . Meyers and Hollowood have performed a detailed 
study of SU(N) gauge theories on S 1 x S 3 with periodic adjoint fermions |119| . In this 
geometry, R 3 is replaced by S 3 , so there are two length scales introduced by the geometry, 
the radius of the three-sphere R = Rs3 and L = R s i. We require min [-R51, Rga] C A so that 
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we are in the weak-coupling region. The projection onto gauge-invariant states, manifested as 
integration over the eigenvalues of the Polyakov loop, ensures non-trivial behavior. Because 
the spatial volume is finite, there is no actual phase transition for finite N, only a crossover 
as R/L is varied. However, the large- N limit does give a phase transition whose behavior is 
closely approximated even for moderate values of N. 

As in the case of R 3 x S 1 , the information about the phase structure is contained in the 
effective potential. Each field adds a factor like 

±Tr R \og[-D 2 -A] (139) 

to the effective potential, where A is the appropriate Laplacian on S 3 for each kind of field, 
± for bosonic and fermionic fields. With Nf flavors of fermion, the effective potential is 
given by 

oo N f 

S(P) = "{(! - z v (nL/R))Tr Adi (P n ) + J2^fWR,m f R)Tr F (P n )} (140) 

n=l U f=l 

where the first term is due to Haar measure. The second term is from the gauge boson, and 
as expected takes the form of a sum over paths winding around S 1 , with z v given by 

°° a p -2L/R _ o p -3L/R 

zJLIR) =2£>+ 2)e-^>/« = (l - - _l% f ■ ("D 
The final term, due to the fermions, has a similar interpretation, with 

oo 

z f (L/R, mR) = 2j2 W + l) e -V(^V2) 2 +m^V«. (i 42 ) 

Although the fermionic contribution cannot be put into a simple closed form, Hollowood 
and Myers derive the useful form 

(j \ 9-m 2 J-} 3 m R f 00 x 2 + - 

m R = —^-K 2 {Lm) - —^-Ki(Lm) + 4 / dx 2 J 4 sin(Lv / a; 2 - m 2 R 2 /R) 
R/L 2 J mR e nx + 1 

(143) 

In the limit m — > 0, this reduces to 

A e -3L/2R 00 

** l i r > °) = (1 _ e -L,R Y - E 2 ^ + i y 

and in the limit R — > oo with fixed m and L, we recover 



-L{i+l/2)/R 



(144) 



z f {L/R,mR)^^^K 2 {Lm) (145) 
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which is the expression that one obtains by working directly on R 3 x S 1 |48| . 
In the case of adjoint fermions, we have 

oo . N 

S(P) = J2~{ 1 ~ z v( nL /R) + N f z f (nL/R,m f R))J2 C0S ( n (. e i ~ 9 i)) ■ ( 146 ) 

n=l *i=l 

In the large- N limit, the effective potential becomes a functional of the distribution p (9) of 
Polyakov loop eigenvalues, with 

S[p(6)] =N 2 fdB f dS' p{6)p{&) fi ~ nL/R ' mK) cos(n(0 - 9')) (147) 

J J n=l n 

where 

f(L/R, mR) = 1 - z B (L/R) + N f z F (L/R, mR) (148) 
and p is normalized to 

d9p(9) = l. (149) 

This expression for S[p(9)} can be made even simpler by writing p in terms of its Fourier 
components 

p( e ) = to E ^ eme ( 15 °) 

ra=— oo 



where p n = TrpP n /N. Then we may write the effective potential as 

N 2 ^ f(nL/R,mR) i 2 



2 ^ n 

n=l 

It is now obvious that the effective potential will become unstable in the contribution of 
iV-ality n precisely when f(nL/R,mR) = 0. These techniques can also be applied to the 
study of gauge theories at finite temperature and density on S 3 x S 1 [1201- 1122] . Another 
interesting geometry is R 2 x T 2 , where SU(N) gauge theories are in the universality class of 
Z(N) x Z(N) spin models because there are two compact directions |123l 1124] . Dimensional 
reduction leads to two-dimensional models where many powerful techniques are available, 
as discussed in Section PVT. 

G. Higgs theories on R 3 x S 1 



As shown by 't Hooft [90, 91j, there is a fundamental conflict between the Higgs mech- 
anism and confinement. The dual superconductor picture of confinement gives a simple 
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picture of this conflict as being between an electric condensate in the Higgs mechanism and 
a magnetic condensate in a confining model. With analytic control of confinement in gauge 
theory on R 3 x S 1 , an adjoint Higgs field can be introduced to study the interplay of con- 
finement and the Higgs mechanism analytically |125[ F126J. The simplest case is SU(2), with 
a classical Euclidean action given by 



with 2 = T 0. The action has a Z(2)h global symmetry given by — > — 0, in addition to 
Z(2)c center symmetry which transform P to —P. 

The scalar field is not gauge invariant, and cannot serve as an order parameter for the 
breaking of the Z(2)h symmetry associated with when gauge interactions are present. 
This is an old problem, a consequence of Elitzur's theorem [127]. Higgs models with scalar 
fields in the fundamental and adjoint representations behave differently. For Higgs models 
with scalar fields in the fundamental representation, the confined and Higgs phases are 
connected [T28] , in a manner similar to the connection between liquid and gas phases. In 
this case, center symmetry is explicitly broken, and large Wilson loops do not have area-law 
behavior due to screening by the scalars. In the adjoint case, center symmetry is preserved 
by the action, and there is a distinct phase transition between the confined and Higgs 
phases. On R 3 x S 1 , there are three distinct gauge-invariant order parameters associated 
with the Z(2)c x Z(2)h symmetry. The first of these is the trace in the fundamental 
representation of the Polyakov loop P itself, (TrpP). It transforms non-trivially under 
Z(2)c but is invariant under Z(2)h- The second is (Trp [P 2 (x) <p(x)]) which is invariant 
under Z(2)c, but transforms non-trivially under Z(2)h- Finally, there is (Trp [P (x) (p(x)}), 
which transforms non-trivially under both groups. 

Using the three order parameters, one can show that there are four distinct phases: a 
deconfined phase, a confined phase, a Higgs phase, and a mixed confined phase. The mixed 
confined phase occurs where one might expect a phase in which there is both confinement and 
the Higgs mechanism, but the behavior of the order parameters distinguishes the two phases. 
In the mixed confined phase, the Z(2)c x Z(2)h global symmetry breaks spontaneously to a 




(152) 



where the Higgs potential is 




(153) 
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Z(2) subgroup that acts non-trivially on both the scalar field and the Polyakov loop. In the 
mixed confined phase, the role of the scalar field is played a linear combination of the Higgs 
field and in the construction of BPS and KK monopole solutions. In all four phases, 
Wilson loops orthogonal to the compact direction are expected to show area-law behavior 
due to unbroken center symmetry in the non-compact directions. This confining behavior 
can be attributed to a dilute monopole gas in a broad region that includes portions of all 
four phases. 

The supersymmetric analog of this model is the Seiberg-Witten model |129j . which is 
an M = 2 supersymmetric gauge theory with gauge group SU(2). Seiberg and Witten 
found that in this model the addition of an M = 1 mass perturbation leads to confinement 
by magnetic monopoles. Recently, Poppitz and Unsal have examined the behavior of this 
model on R 3 x S 1 , and concluded that the confined phase seen for small compactification 
circumference on R 3 x S 1 is connected to the confining phase at infinite compactification 
circumference [130J. In their work, Euclidean monopoles in which a linear combination of 
A4 and plays the role of the scalar field appear in a very similar fashion to the non- 
supersymmetric model. 

V. CONFORM ALITY AND DUALITY 

A. Non-trivial fixed points in d = 3 + 1 gauge theories 

Up until this point, we have implicitly restricted ourselves to theories like QCD which 
have a simple renormalization group structure: they have an ultraviolet fixed point at g — 0, 
and hence are asymptotically free, and an infrared fixed point at g = 00 when the theory is 
defined on R 4 . Gauge theories with non-trivial fixed point structure are known to occur if 
the number of particle representations included in the theory is sufficiently large. The first 
two coefficients of the perturbative contribution to the renormalization-group /3 function of 
a gauge theory are independent of the renormalization scheme, and given by 



dg 2 





(154) 



d log q 2 
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where 



bx = jC 2 (G)-^N f T(R) 
34 

b2 = y [C 2 {G)f ~ N f T(R) 



yC 2 (G)+4C 2 (i?) 



The coefficients b\ and 62 depend on the group and the group representations of the particles 
via C2(R), the value of the quadratic Casimir operator in representation R {G denotes 
the adjoint representation, so C^G) = N ), while T(R) is the normalization of the group 
generators in R: Tr {T%T%) = T(R)5 ab . 

Three possible behaviors may be obtained from perturbation theory: 

1) bi < : the theory is no longer asymptotically free at high energies, and g = is 
an infrared fixed point, as is the case in QED. For a non-Abelian gauge theory, 
this is typically the case when the matter content of the theory becomes too 
large. Only the gauge fields themselves give a negative contribution to b\, so a 
sufficient number of additional matter fields can drive hi negative. 

2) b% > and 62 > 0: in this case the coupling runs out of the region of pertur- 
bative computability in the IR. This is the QCD-like case, where the theory is 
asymptotically free at high energies (the UV) and confining at low energies (the 
IR). 

3) bi > and 6 2 < 0: Perturbation theory predicts a non-trivial fixed point, given 

by 

g 2 BZ = -I67A (155) 



These three behaviors are shown in 10 The corresponding phase structure as a function of 



the number of flavors is shown in figure [TTJ As was first pointed out by Banks and Zaks 
[131| . this fixed point may or may not be located in the region where perturbation theory 
is reliable. However, when it is located in that region, the prediction of an infrared-stable 
fixed point is reliable. Note however that a change in the renormalization prescription, which 
redefines the coupling constant, will change the beta function. This in turn will change the 
location of a non-trivial fixed point. Given that such theories are possible, we would like 
to understand their properties. Such theories are not generally candidates for describing 
the real world. Theories with non-trivial infrared fixed points have correlation functions 
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Figure 10: Three possible behaviors for the beta function. If Nf < Nf, the theory is asymptotically 
free but has no non-trivial fixed point. If Nf > N^ F , the gauge theory is no longer asymptotically 
free. For N* f < N f < Nf F , there is a conformal window with a non-trivial infrared fixed point. 

that decay algebraically at large distances, in a manner familiar from critical behavior at 
second-order phase transitions. Requiring that models not have an infrared fixed point is 
a condition that likely requires non-perturbative understanding that lattice gauge theories 
can supply. The case of a non-trivial infrared fixed point, case 3), is intermediate between 
the familiar case 2) and non-asymptotically free theories. Typically, for a given gauge group 
and choice of additional representation R, there is some number of flavors Nf- F such that 
asymptotic freedom is lost, and some smaller value Nf where an infrared fixed point appears. 
Theories with N* f < N f < Nf F are said to lie in the conformal window. 

There is a close relation between the existence of the conformal window and possible 
beyond the standard model physics associate with the Higgs field. In the standard model, 
the Higgs field is a fundamental scalar field, in the sense of not being a composite; in 
fact it is the only fundamental scalar field in the standard model. Because the bare Higgs 
mass receives corrections which are quadratic in any ultraviolet cutoff, there is a problem 
of "naturalness" or "fine-tuning" in achieving a sufficiently small Higgs mass. This problem 
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Figure 11: Phase structure of a theory with a conformal window as a function of the number of 
flavors Nf. 

would be circumvented if the Higgs were a bound state of some new fermions rather than 
a fundamental scalar. Theories of this type that have a non-fundamental Higgs are called 
technicolor models, and the gauge group of the new fermions is the technicolor group; see 
|132j for a recent review. In the standard model, the Higgs is responsible not just for 
the gauge boson masses, but also for the fermion masses. The class of models needed for 
successful phenomenology are called extended technicolor models, and seem to require a 
gauge theory very close to the conformal window. Such theories are often referred to as 
"walking" because the running coupling constant moves slowly if the beta function is small, 
as it would be near a zero. 

There has been much work done on the location of the conformal window but definitive 
results have not yet been achieved. Much of the analytic work uses the chiral order parameter 
ijjip as the order parameter for the transition from a confining phase, where ipip is expected to 
be non-zero, to a conformal phase which has no dynamically generated scale. One method 
is based on truncated Schwinger-Dyson equations |133H138] while another approach uses 
renormalization group methods |139| . In many supersymmetric theories, the boundary of 
the conformal window can be determined, as discussed in |140J. There are also estimates 
for non-supersymmetric models based on ideas from supersymmetry |141j . Introductory 
reviews of lattice gauge theory research in this area can be found in |142H145] . Recent 
research includes the study of model with fermions in the fundamental representation |146F 
150J, the adjoint representation |151H157] , and the symmetric representation |158f Tl62j . 

B. Conformality and duality in Abelian models 

There is an interesting class of lattice models that share some of the features of gauge 
theories with conformal windows: Z(N) spin systems in two dimensions and Z(N) gauge 
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theories in four dimensions. Our understanding of these models originates in the two- 
dimensional XY model, which is a spin model with a global U(l) symmetry. The order 
parameters of the model are the operators 



where (x) takes on values between and 27r. The XY model provides the principle example 
of a phase transition driven by topological excitations, the Berezinsky-Kosterlitz-Thouless 
(BKT) transition j!63l 1164] . Low-temperature arguments indicate that the XY model has 
a low-temperature gapless phase represented by a line of critical points, along which critical 
indices vary continuously with the temperature. On the other hand, high-temperature ex- 
pansions indicate a gapped phase at high temperatures. The phase transition that separates 
the two phases is driven by vortices, configurations of spins that have a non-trivial winding 
number. In a continuum notation, these are configurations that have 



where n is a non-zero integer. Vortices act as a two-dimensional classical Coulomb gas, 
interacting via a long-ranged logarithmic interaction. In the low-temperature phase, vortex 
pairs bind tightly to form bound states of zero vorticity and have no effect on the large- 
distance behavior. At the critical temperature, vortices unbind to form a classical Coulomb 
plasma, giving rise to a mass gap, the Debye mass of the plasma. From a naive continuum 
field theory point of view, the XY model has only massless spin-wave excitation, with an 
Euclidean action given by 



where J 2 is a coupling constant inversely proportional to the temperature. It is convenient 
to rescale the field by — > <p/J so that the kinetic term of the action has conventional 
normalization. Vortices are created and destroyed by operators of the form 




(156) 




(157) 




(158) 




(159) 



where <p and <p are related by 



(160) 



The effect of vortices gives rise to an effective action 




(161) 



67 



yi 



irrelevant 



t t t 



relevant 

AAA 



-► j 



-2 



Massless n /2 Massive 
spin-wave BKT vortex 
phase transition phase 

Figure 12: Relevance of V\ for the sine-Gordon model. 

where A is a cutoff on the order of the lattice spacing and y p is a dimensionless activity for 
a vortex of winding number ±p |165l 1166] . This generalized sine-Gordon model gives rise to 
Coulomb gas of vortices when the partition function is expanded in the activities. Generally 
speaking, it is only necessary to consider the p — ±1 vortices. This model can be perturbed 
to a model with Z(p) symmetry by the addition of an operator — h p A 2 cos which 
explicitly breaks the U(l) symmetry down to Z(p): 

~J 2 



S, 



eff 



d x 



{d<f>y - 2 Vl A 2 cos 2tt J0 - 2h p cos (p<j)/J) 



(162) 



This model has a duality under the interchange 

<H- (f) 

h p -H- yi 
P 



J 



O 2nJ 



which is a generalization of Kramers- Wannier duality for the Ising model |167[ 168J. 

A simple understanding of the renormalization-group behavior of the XY model is pro- 
vided by analyzing the multiplicative renormalization of the spin-wave and vortex operators 
associated with normal ordering |169[ 1170] . For a massless free field, the relation between 
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Figure 13: Actual sine-Gordon flow. 



the bare field exp (i/3<p) and the normal-ordered field N [exp (i/3(j))] is 



2Jt 



9 3 



[e iW ] 



(163) 



where /i is arbitrary. An intereaction term of this type will be relevant if /3 2 < 8tt and 
irrelevant if /3 2 > 8tt. When applied to the XY model, this implies that vortices are relevant 
only if J~ 2 > 7r/2, in the high-temperature phase. This leads to the simple picture shown 
in Fig. [12} Taking into account the renormalization of the vortex activity y± leads to a 
more complete picture of the renormalization group flow for the XY model, as shown in 



Fig. 13 Applying the same simple condition to the Z(p) model indicates that the spin- 
wave symmetry breaking term S p + S_ p is relevant if J~ 2 < 8ir/p 2 , corresponding to low 
temperatures. If p > 4, there is a gap between the regions where neither spin-waves induced 



by hp 7^ nor vortices are relevant, as shown in Fig. 14 In this intermediate conformal 



window, correlation functions decay algebraically; there is no mass gap. This phase structure 
is only part of a larger picture for more general Z(p) lattice models. As we have seen in 



Sections III and |I V| systems with Z(p) symmetries can break spontaneously to a non-trivial 



subgroup of Z(p). This also occurs in d = 2 Z(p) spin systems |171H173] ; for example Z(6) 
can break spontaneously to Z(2) or Z(3). For p = 2 and p = 3, the regions of vortex and 
spin-wave relevance overlap, and must be handled as special cases. 

These results extend to gauge theories in several different ways. The analysis of two 
dimensional systems by Jose et al. [165J can be extended to Abelian lattice gauge the- 
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Figure 14: Relevance of and V\ operators in a model for p > 4. 

ories. In four-dimensional Abelian lattice models, the analogs of the point-like vortices 
of the two-dimensional spin models are closed magnetic monopole loops j!74J . Similar to 
two-dimensional spin systems, four-dimensional Z(p) lattice gauge theories have a duality 
between monopole loops and loops of charged particles. For p > 4, there is a gapless in- 
termediate phase where neither the electric nor the magnetic operators are relevant |175| . 
In three-dimensional Abelian lattice models, the topological excitations are point-like, as is 
the case in Polyakov's treatment of the continuum three-dimensional SU(2) gauge theory 
broken to U(l). This is not accidental, because there is a clear relation between Abelian 
duality on the lattice and in the continuum. At finite temperature, i.e., in a R 2 x S 1 ge- 
ometry, Polyakov's model can be analyzed in a manner analogous to d = 2 spin systems. 
Three-dimensional monopole solutions made periodic on S 1 play the role of vortices, and 
early work identified the phase transition to be in the XY-model, or U(l), universality class 
|176| . As in 3 + 1 dimensions, adoint Polyakov lines associated with W ± gauge bosons ensure 
that the symmetry group of the low-energy theory remains Z(2) [177]. 

Another interesting connection between the behavior of d = 2 spin models and d = 4 
gauge theories is found in the suggestion by Kaplan et al. that there may be two nontrivial 
fixed points in gauge theories with a conformal window |178| . In this scenario, there is a UV 
fixed point at g 2 = 0, followed by a non-trivial infrared fixed point, followed by a non-trivial 
ultraviolet fixed point at a larger value of the g 2 . As Nf decreases, the two non-trivial fixed 
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points annihilate, which is to say they move off into the complex plane. This behavior would 
give rise to BKT scaling behavior in the region near where the fixed point merger occurs. 



C. Relevance of topological objects on R 3 x S 1 

Unsal and Poppitz have proposed an alternative criterion for locating the boundary be- 
tween confining and conformal gauge theories |179[ 1180] . Their work is based on the relevance 



or irrelevance of the topological objects relevant for confinement on R 3 x S 1 , analogous to 
the criterion for the relevance of vortices in d = 2 Abelian spin systems. For all center- 
symmetric gauge theories on R 3 x S , the mass scale induced by topological excitations can 
be written in the form 

1 r s^ 2 i 

(164) 



m (L) ~ — exp 
ij 



-Q- 



y (l)n\ 

where q is O (1) and depends on the details of the theory. With the use of the renormalization 
group, we can write this as 

m{L) ~ K(KL) ql1 - 1 (165) 

where b\ = b\/N / Requiring that the theory remain in the center- symmetric, confined phase 
as L becomes large puts a constraint on the theory. In the limit L — > oo, the topological 
excitations are relevant if qb\ > 1 and the theory is confined for large L. If qb\ < 1 , the 
theory is not confined. While conceptually simple, the details of the calculations for different 
models are somewhat technical. The parameter q is an integer determined essentially by 
how many monopoles, each of action 87r 2 /g 2 N, are required in the process that generates 
the mass gap. For flavors of Dirac fermions in the fundamental representation, the 
conformal window is estimated to lie in the range 

< Nf < y n ( 166 ) 
while for Nf flavors of Weyl fermions in the adjoint representation the range is 



4 < Nf < y (167) 



where = 2Nj° if Dirac fermions are used. Estimates of the conformal window for 

other representations can also be made [179J. The estimates are comparable to estimates 
made using other methods, but the connection between this method and others is at present 
unclear. 
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VI. GAUGE THEORIES IN THE LARGE-N LIMIT 



Many theories simplify in the large- iV limit, where N is the dimension of some internal 
symmetry group such as SU(N) taken to be large. Although physical QCD has N = 3, there 
is reason to believe the large- N limit provides a good description of many features |181[I182| . 
Lattice simulations at finite values of N have been able to provide accurate extrapolations to 
the large- N limit for some quantities |183[ 1184] . For example, lattice simulations of SU(N) 
gauge theories in (2 + 1) dimensions have provided evidence for 1/N corrections in string 
tension scaling laws |3H]- The AdS/CFT correspondence, which relates gauge theories to 
gravity duals in the large- N limit |185] , provides another reason for interest in large- N gauge 
theories; see |186l 1187] for reviews. 

The large- N limit began as an approximation in statistical mechanics [188] . but was 
quickly applied to models with scalars cf) a [189] or fermions ip a |190] in the vector represen- 
tation of the groups O (N) or SU (N) . Such models are exactly solvable in the large- N 
limit. Field theories with particles in the adjoint representation, such as gauge bosons, are 
much more difficult. 

The original approach of 't Hooft to the large- N limit in SU(N) Yang-Mills theories was 
based around a detailed analysis of classes of Feynman diagrams using an ingenious double 
line notation |181j . This notation is based on the group properties of SU(N) propagators: 
for fundamental representation fermions, we have 

(r(x)A(y)) = S a b S(x-y) (168) 

but for gauge bosons, which are in the adjoint representation, we have 

1 



A« b (x)At d (y)) 



A a A c A a A c 

°d°b ]\[ b d 



D^ix-y) (169) 



In this notation, quarks in N, the fundamental representation of SU(N), and antiquarks in 
N are represented by single lines with arrows. Gauge bosons in the adjoint representation 
are represented by a double line with one arrow in each direction; other representations may 
be represented by other combinations of arrows. 

The Euclidean Lagrangian for an SU(N) gauge theory can be written as 

Le = ^ 2 Tr F [F%\ = ^Tr F [i^J (170) 
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we have introduced A = g 2 N, which will be taken to be finite in a limit where g — > as 
iV — > oo. It is obvious that every gauge boson propagator will carry a factor of 1/N, and 
the cubic and quartic vertices will carry a factor of N. Fermion fields can be rescaled by 
ip — > N 1 / 2 ^ so that similar properties hold. 

Any vacuum diagram can be regarded as a polygon, where every closed loop can be 
regarded as a face. A vacuum diagram can then be characterized by its number of faces 
F, number of vertices V and number of edges E. Each face generates a factor of iV from 
a trace. Every edge is a propagator, and carries with a factor of iV -1 while every vertex 
carries a factor of N. Thus each vacuum graph carries a factor of 

N F-E+V (m) 

However, F — E + V is the Euler number xe of the surface, a topological invariant. It may 
be calculated from the number of handles if of a surface and the number of holes B (for 
boundary) as 

X e = 2-2H-B. (172) 

Thus the maximum power of N associated with a vacuum diagram is N 2 , for diagrams with 
H = B = 0. Because each fundamental representation loop gives rise to a hole, fundamental 
representation particles begin contributing at order N. This behavior is what we see in the 
perturbative calculation of the effective potential for the Polyakov effective potential. The 
effects of quarks are suppressed by a factor Nf/N relative to gluons. It should be noted 
that there is another interesting large- N limit, the Veneziano limit, in which Nf is taken 
to infinity at the same time as N, with the ratio x = Nf/N fixed [191J . This generates a 
large- N expansion where particles in the fundamental representation are not suppressed. 

A key feature is large- iV factorization. Suppose, for example, that 0\ and O2 are operators 
transforming according to the adjoint representation of the group. Then TrOi/N has a 
smooth limit as N — > 00. Large- N factorization states that in the large- N limit 

(^ Tr0 4 Tr ° 2 ) -+ ( A7 Tr01 ) ( A7 Tr ° 2 ) + ° ( 1/iV2 ) (173) 
up to terms of 0(1/N 2 ). For Wilson loops, large- iV factorization leads to a closed, exact 
equation in the iV — > 00 limit, the Migdal-Makeenko equation j!92j . It was quickly realized 
that many aspects are better captured by postulating the existence of a large- iV saddle point 
of some large- N effective action analogous to the gap equation in vector-like models. This 
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saddle point, really a gauge orbit in the space of fields, is often referred to as the master 
field |193| . While it is possible to find the master field for some model systems in lower 
dimensions with fields in the adjoint representation [194], the problem of determining the 
master field for four-dimensional gauge theories has not been solved. 

In recent years, much of the work on gauge theories in the large-iV limit has em- 
ployed gauge-string duality to explore strong-coupling behavior using generalizations of the 
AdS/CFT correspondence [185J. For a review emphasizing the relation of gauge/string du- 
ality to finite temperature QCD, see P-87J. Attempts to model the phase structure of QCD 
based on the AdS/QCD approach [J95] may be found in [196- 1199] . 

A. Eguchi-Kawai models 

A direct attempt at constructing a large- iV reduction for lattice gauge theories was first 
made in 1982 by Eguchi and Kawai pOOj. They defined a lattice gauge theory on a single 
lattice site and showed that Wilson loops in an infinite volume theory could also be obtained 
from their reduced model in the large- A limit if two conditions were met: large- A factor- 
ization and center symmetry must both hold in all directions. This latter restriction implies 
that the reduction can only hold in the confined phase. 

The original Eguchi-Kawai (EK) model may be derived by observing that the master 
field must be translation invariant. For a lattice gauge theory, this means that the value of 
the link variable at every site must be a gauge transform of the field at the origin 

U IM (x)^D\x)U fM D(x + fi) (174) 

where denotes the link variables at the origin. The matrix D(x) can be taken to be a 
diagonal matrix of the form 

D(x) = exp {-ix ■ P) (175) 

where P M is a diagonal Hermitian matrix of the form P^ = diag (pi^, ■■■,Pnh)- Gauge- 
invariant quantities immediately collapse. For example, a rectangular RxT Wilson loop in 
the \x — v plane becomes 

(Tr (U,) T (U„)\UI) T (U!) R ^ (176) 
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and the Wilson action collapses to 



Sek = ^ ReTr [u&Uffi] (177) 

where VT is the volume of space-time. Note that global center symmetry transformations 
of the form — > z^U^ is a symmetry of Sek- 

In the infinite volume gauge theory, the expectation value of non-closed paths such as 

Tr(u^ T \ (178) 

is zero as a consequence of gauge-invariance and Elitzur's theorem. In the Eguchi-Kawai 
model, this expectation value will vanish if center symmetry is unbroken, but need not vanish 
if center symmetry is unbroken. A detailed comparison of the Schwinger-Dyson equations 
for both theories confirms that they are equivalent only if center symmetry is unbroken. The 
similarity with the case of finite temperature immediately suggests that center symmetry will 
be broken in the large- N limit, and the equivalence will fail. This turns out to be the case 
|201H203] . Several attempts have been made to modify the original EK model to maintain 
center symmetry. The first was the Quenched EK (QEK) model, in which the eigenvalues 
of the link variables are treated as quenched variables, forcing the link variables to maintain 
center symmetry |201] . While early lattice simulation did indicate that the QEK model 
was viable, but recent lattice simulations provides strong evidence for the breakdown of the 
reduction |204J . Another variant of the original model that originally seemed promising is 
the Twisted EK (TEK) model |205| . but more recent extensive lattice simulations indicate 
center symmetry breaking in this model as well |2U6H2U8] . An interesting alternative to 
single-site reduction is to simulate large-iV theories on lattices of size N*, and keep N s 
sufficiently large that the deconfinement transition is avoided |184[ [209J . As long as center 



symmetry is maintained, such models will be equivalent to infinite-volume theories in the 
limit of infinite N. This behavior, a hallmark of a successful large-iV reduction, is often 
generally described as volume independence. 



B. Planar equivalence 

The most prominent example of planar equivalence is the equivalence, at large N, of 
the M = 1 super- Yang-Mills model to an SU(N) gauge theory with fermions in either the 
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symmetric or the antisymmetric two-index representation. In other words, in the large-A 
limit, SU (N) gauge theories with fermions in the adjoint representation (Adj) are equiva- 
lent to the same gauge theory with fermions in either the symmetric (S) or antisymmetric 
(A) representation. The history of the subject has an interesting arc; see |210] for an intro- 
ductory discussion. There is a very general formalism for describing many forms of large- iV 
equivalences |211j . with which a detailed discussion of many forms of planar equivalence 
can be carried out |212H214] . For our purposes, the equivalence of SU(N)Adj t A,s can be 
understood from the equivalence of the loop equations for the different theories |215| . Su- 
persymmetric gauge theories have Majorana fermions in the adjoint representation which 
are superpartners of the gauge bosons. While the effects of fundamental representation 
fermions are suppressed in the large-A limit, the effects of adjoint representation fermions 
are not, because the number of components in the representation grows as N 2 in the large- 
A limit. In general, the dimensionality da of a given representation R is easily determined 
from its character Xr(U) = TtrU via cLr = Xr(-0- F° r example, we have for the adjoint 
representation 

Xa#(U)= X f{U)xp(U)-1 (179) 

so d,Adj — A 2 — 1. The characters of the symmetric and antisymmetric representations 
Xs/a{U) are formed from the product of two characters in the fundamental representation 
and given by 

Xs/a(U)= 1 -[ X 2 f (U)±Xf(U)] (180) 

with dimensions N(N ± l)/2. These representations are sometimes referred to as bifermion 
representations. Because the number of fields grows as A 2 , we expect the effects of internal 
loops of S/A fermions to survive in the large- N limit in the same way that adjoint fermion 
loops do. 

Let us imagine integrating out the fermions in an SU(N) gauge theory with either adjoint, 
symmetric or antisymmetric fermions [216]. In all cases, the logarithm of the fermion deter- 
minant can be expressed as a sum of Wilson loops in R 4 ; in a geometry like R 3 x S* 1 or T 4 , 
the sum includes Polyakov loops as well as Wilson loops. In the large- N limit, factorization 
tells us that on R 4 we have 

2 (xs/a(W)) ~ ( X f{W)) ( X f(W)) ~ ( X f(W)) (x* f (W)) ~ ( X A dj (W)) (181) 
where the expectation value is over the gauge fields. The apparent mis-match of a factor 
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of two disappears in the functional determinants because the fermions in the symmetric or 
antisymmetric representations are Dirac fermions, while the adjoint fermions are Majorana 
fermions with half the degrees of freedom. A proof of the equivalence at the level of Feynman 
diagrams can also be given using 't Hooft's double line notation [217J. In the confined phase 
on T 4 or a similar geometry, all the expectation value for the Polyakov loops are zero 

2 ( X s/a(P)) ~ (xf(P)) (Xf(P)) ~ (Xf(P)} (x* f (P)) ~ (XA dj (P)} = (182) 

and any of the reduced gauge theories that is confining is completely equivalent to the 
theories on R A , realizing volume independence. However, this equivalence fails in non- 
confining phases, because the iV-ality of the symmetric and antisymmetric representations 
is two, but for the adjoint representation it is zero. Thus we return once again to the issue 
of maintaining confinement in small volumes. 



Confinement for finite volume 



As we have seen in section IV, suitable modifications of the action can lead to a confining 
phase on R 3 x S l for small circumference L. It is natural to look at these methods for the 
possible construction of useful reduced large-iV models. However, the confining phase is 
only one of many possible phases as N becomes large |4"8l 149] . In principle, it is possible to 
modify the large- N gauge action on T 4 by adding double-trace deformations 

[N/2] 

S^S+ J2J2 ak / d 4 xTr A P tl {x) k . (183) 

k=l p. J 

where (x) is the Polyakov loop in the direction fi based at the spacetime point x. There 
is a harmless redundantcy induced by the integration over x, because all points a given loop 
path are integrated over. For a single-site lattice action, the corresponding modification is 

[N/2] 

s^s+j2J2 akTrAU t- ( 184 ) 

k=l fi 

because the single-site link variable is the Polyakov loop. If the coefficients are sufficiently 
large, then the system will remain in the confined phase, with (TrAP k ) = . In the confined 
phase, the precise values of the coefficients don't matter; we would have a self-consistent 
finite- volume reduction of the pure gauge theory |93j. 
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Adjoint fermions with periodic boundary conditions offer another possible route to the 
construction of a large- N finite- volume reduction [971 ESI 11191 1218H226] . This model is also 
of direct phenomenological interest. In physical QCD with iV = 3, we normally think of 
quarks in the fundamental representation as playing no role in the large- N limit. However, 
antiquarks in physical QCD are in the 3 representation, which is an antisymmetric represen- 
tation. In this sense, physical QCD has two different large-iV limits. The large- N limit, an 
SU(N) gauge theory with Nf Dirac fermions in the adjoint representation is equivalent to 
the theory with 2Nf Dirac fermions in the antisymmetric representation, which is QCD with 
2Nf flavors when N = 3. Analytic calculations of the lattice form of the one-loop effective 
potential [971 11191 12181 12211 1222] indicate a rich phase structure as the fermion mass is varied, 
similar to that found in the continuum on R 3 x S 1 |18| I119j . However, it should be kept in 
mind that chiral symmetry breaking effects, which affect Polyakov loop terms in the effective 
potential, are generally not included in these calculations; at present, these effects can only 
be estimated using phenomenological models [71J . Thus lattice simulations are essential not 
only to determine the properties of the confined phase, but also to confirm that it exists. 
The key issue is finding a confined phase in lattice simulations that survives the extrapo- 
lation to the continuum limit where the bare lattice coupling g 2 goes to zero. Of course, 
the success of the large- N reduction may depend on the lattice fermion implementation. In 
the single-site model, analytic calculations indicate that center symmetry is spontaneously 
broken with naive fermions but unbroken with overlap fermions |223j . Lattice simulations 
of SU(3) with 2 flavors of staggered fermions in the adjoint representation on R 3 x S 1 show 
that a confined phase does exist for sufficiently small fermion mass [96J. Simulations of the 
single-site theory using Wilson fermions also show a region in the parameter space of fermion 
mass and lattice coupling where center symmetry is unbroken and an extrapolation to the 
continuum limit appears plausible [98, 225j;see also |224j . Simulations with overlap fermons 
[226J also indicated a confined region. The most recent lattice results [227J of a single-site 
model with two flavors of Wilson fermions show a center-symmetric region for N as large 
as 53 and 't Hooft coupling g 2 N as small at 0.005, so the prospects for a successful large- N 
reduction appear bright. 
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VII. CONCLUSIONS 



We have entered a new phase in our understanding of the phase structure of gauge 
theories. The phase structure of gauge theories at finite temperature was known to be rich. 
By itself, however, temperature alone does not give us a large window where continuum 
analytical results and lattice simulations are both useful. We have now a class of theories 
that can be successfully studied using both analytical methods and lattice simulations, built 
on our understanding of finite temperature physics. On _R 3 x S 1 , the use of double-trace 
deformations or periodic adjoint fermions allows us to study confinement in a region where 
semiclassical methods are valid, and check our results with lattice simulations. At the same 
time, we have access to many new phases, most of which are partially-confining phases. 
Nor are we restricted to an R 3 x S* 1 geometry; other geometries are available, and largely 
unexplored. By extending finite-temperature physics to this larger class of models, we have 
seen how instantons and non-Abelian monopoles can play a role in confinement, realizing 
some long-held ideas about the nature of the confining phase. At the same time, we have 
seen how topological effects can survive in the large- N limit. We also have a new, promising 
class of single-site large- N models to explore. 

There is more to be done. Most of the lattice simulations to date have only checked 
analytic predictions of overall phase structure, and predictions for, e.g., the behavior of 
string tensions have not been checked in simulations. Although semiclassical techniques 
have a natural range of validity on R 3 x S 1 given by NAL < 1, we do not know the 
general validity of the concepts useful for small L. We know that a double-trace deformation 
connects the small- L confining phase to the large- L confining phase without an intervening 
phase transition, but we do not know how physical quantities change along a path that 
connects the two. We do not yet have a complete understanding of string-tension scaling 
laws, even in the small- L region, and there are many other questions that remain to be asked 
and answered. 
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